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Abstract 


The  purpose  of  this  work  was  to  develop  a dynamic  partial  order  model  of 

It 

electrolyte  solutions  which  avoids  inconsistencies  present  in  the  Debye-Huckel 
and  in  related  theories.  It  is  shown  that  our  results  are  consistent  with  activity 
coefficient  data  in  aqueous  electrolyte  solutions.  The  main  features  of  our  work 
are  a partial  long-range  order  which  increases  with  concentration  due  to  the 
increase  in  the  coulombic  energy  of  interaction  and  decreases  with  temperature, 
a contribution  from  ionic  fluctuations  because  it  is  shown  that  displacements 
from  average  positions  lead  to  a symmetric  increase  in  energy  which  is  not 
cancelled  out,  and  a hydration  correction  for  the  effect  of  hydration  upon  the 
effective  ionic  concentrations.  Ions  are  given  partial  counter- ion  character 
as  a function  of  the  degree  of  order  so  that  the  free  energy  of  ionic  interactions 
can  be  calculated  rigorously,  as  if  the  lattice  were  perfect,  without  resort  to 
cluster  integrals  in  the  case  of  single  salt  solutions. 
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Introduction 

Our  purpose  in  developing  this  partial  order  theory  of  activity  coefficients 
in  electrolyte  solutions  was  to  bypass  inconsistencies  present  in  the  Debye- 

ii 

Huckel  theory  and  in  related  ones  such  as  extensions  applied  to  specific  ionic 
interactions.  These  inconsistencies,  some  of  which  have  not  been  pointed  out 
before,  arise  primarily  from  the  use  of  the  Boltzman  distribution  to  a problem 
the  very  essence  of  which  is  the  interaction  of  charged  particles,  and  from 
neglecting  the  quasi-lattice  nature  of  ionic  solutions  which  precludes  net 
electrical  densities  from  decaying  monotonically  with  distance  from  a central 
(reference)  ion. 

We  will  start  this  work  with  a brief  sketch  of  the  derivation  of  the  Debye- 

II 

Huckel  equation,  followed  by  a critique  which  will  bring  to  light  logical  impasses 
not  previously  revealed. 

We  will  then  examine  early  lattice  models  in  which  a fixed  and,  generally, 
perfect  degree  of  long  range  order  was  assumed  causing  the  resulting  equations 
for  activity  coefficients  to  only  apply  over  limited  ranges  of  concentrations. 

Finally  we  will  present  our  model,  valid  at  all  concentrations,  in  which 
the  degree  of  order  increases  with  concentration  because  of  the  increase  in  the 
coulombic  energy  of  interaction  and  decreases  with  increasing  temperature. 

Our  model  will  be  developed  step-by-step  for  greater  clarity  with  the  successive 
addition  of  relevant  factors,  namely,  partial  long-range  order,  translational 
energy  and  hydration.  The  translational  energy  is  especially  important  because 
oscillations  of  ions  around  central  positions  lead  to  a symmetric  energy  contri- 
bution which  is  not  canceled  out. 
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We  will  not  Drove  conclusively  that  our  model  is  the  only  right  one 
because,  as  in  all  other  models  of  activity  coefficients,  we  have  an  adjustable 
parameter  of  suspected  physical  significance  but  of  unknown  a priori  analytical 
form  as  a function  of  the  concentration,  but  we  will  show  that  the  model  is 
plausible  and  is  consistent  with  experimental  data.  Thus,  our  main  purpose 
at  this  stage  is  explanatory  rather  than  predictive,  a feature  common  to  all 
models  of  activity  coefficients  which  require  empirical  fitting  at  some  stage  in 
their  application. 
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1 . Debye-Huckel  Theory  and  Critique 
It  will  ease  our  critical  review  and  refresh  the  memory  of  some  of  our 

ii 

readers  if  we  sketch  the  derivation  of  the  Debye-Huckel  (1923)  equation  for 
activity  coefficients  of  electrolytes  in  solution  at  the  outset  of  this  work. 

ii 

a)  The  Debye-Huckel  equation 

Consider  a central  (reference)  ion  J which  is  assumed  to  sense  the  other 
ions  present  in  the  solution  as  a smeared  continuous  distribution  of  charges 
with  a charge  density  p in  a volume  element  dV  placed  at  a distance  r 

( 3 ) ( J ) 

from  J.  Poisson's  equation  at  dV  is  then 


2„  4tt 

7 0(J)  = 'D  C(J) 
e 


(j- j) 


in  which  0 is  the  electrical  potential  at  dV  resulting  from  the  fields  produced 
( J ) 

by  all  the  ions  in  the  solution  and  D is  the  dielectric  constant. 

e 

Assume  next  that  the  Boltzman  law  is  valid  for  the  distribution  of  ions 
around  J and  let  i represent  each  icnic  species  present  in  the  solution  (e.  g.  , 
i = Na+,  Cl  in  a NaCl  solution).  Then,  one  may  write 


/ Z-e0/T\ 

> ( J ) , 

n.  = n.  exp  (-— — “) 
1 1 kBr 


(1-2) 


n.  is  the  number  of  ions  of  type  i per  unit  volume  at  dV,  n.  is  the  bulk  value  at 

great  distances  from  J.  z.  is  the  valence  of  i,  c is  the  electronic  charge,  k„  is 
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Boltzman's  constant,  and  T is  the  absolute  temperature. 
The  net  charge  density  at  dV  is 


(J) 


= I a,  z.e  - £ n.z  e exp  (- 

• i i ‘ii  1 ' 


Z €^/  T \ 

_} ill 


1 1 


k T 
B 


(1-3) 


which  decays  to  the  bulk  value 


ti) " ^ Vie  “ 0 


(1-4) 


at  large  distances  from  J (and  implicitly  from  other  ions  so  that  0 may  be 

( J ) 


essentially  zero). 

The  intrinsic  concept  at  this  point  is  that,  because  of  coulombic  attraction, 
there  is  an  excess  of  counter-ions  around  J,  i.  e.  , of  ions  of  opposite  charge, 
and  a defficiency  of  ions  of  like  charge  but  that  one  reaches  electrical  neutrality 
when  the  volume  element  dV  is  moved  far  enough  from  J. 

By  introducing  Eqn.  (1-4)  into  (1-1)  Poisson's  equation  becomes 


2„  4tt  ~i’"m 

7 0U)  = -5~  f niV  “■>  (-iT7“ 


z.e0, 


d-5] 


B 


There  are  two  difficulties  with  this  equation;  its  rigorous  solution  is  to  be 
avoided  and  it  is  not  compatible  with  the  principle  of  the  superposition  of  fields 
In  effect,  if  the  ionic  charges  were  doubled  one  should  expect  the  potential  at 


1 -3 

dV  to  also  be  doubled  but  this  does  not  occur  due  to  the  exponential  nature  of 
Eqn.  ( 1-3). 

tf 

Debye  and  Iluckel  avoided  these  problems  by  expanding  the  exponential 
term  by  means  of  the  series 


and  obtained 


(1-6) 


n z e 
i i 


[1 


!if!m  +_l 

kBT 


Z G0/  T, 

■ i-JJ)) 

kBT 


(1-7) 


The  first  term  in  the  expansion  is  zero  because  of  electrical  neutrality  (see 
Eqn.  (1-4))  and  it  can  be  shown  that  all  the  terms  in  even  powers  vanish  for 
symmetric  (but  not  for  asymmetric) electrolytes. 

Eqn.  (1-7)  is  then  reduced  to 


where 


2„  4rre  _ 2„  2„ 

V 0/ti  = n l x En.z.  0 = K 0 

(J)  D k T ill  (J)  (J) 

e Jo 


(1-8) 


2 2 
_ 4ne  En.z. 

V2 L-  1 1 

K D k T 

e D 


8TTe2l^c>N 
S A 

1000  Dk,T 
e B 


I 


(1-9) 
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( c) 

1,  being  the  molar  ionic  strength  of  the  solution  and  N Avogadro's  number, 
S A 

Converting  Eqn.  ()-8)  into  polar  coordinates  and  integrating,  one  obtains  the 
solution 


-Kr 


(J) 


K r 


(J) 


% ’ ki 


+ k 


(J) 


2 r 


0-10) 


(J) 


When  r tends  to  infinity,  0 must  remain  finite  and,  therefore,  k =0. 

( J ) ( J)  2 

Electrical  neutrality  requires  that 


I 


4',0(j)dr  = -v 

D 


0-11) 


where  a^  is  the  distance  of  closest  approach  of  the  center  of  other  ions  to  the 
center  of  J.  k^  is  calculated  by  introducing 


■Kr 


(J) 


0 = k 

v/  Tt 


(J)  1 r 


0-12) 


(J) 


into  the  expression  (1-3)  for  the  charge  density  after  it  is  linear ized, and  inte- 


grating by  parts.  One  then  obtains 


*ar^  -*r, 

z e D (. 

P,J)  = 1T  ^ iT-,aDx£T^ 


(1-13) 


K can  be  shown  to  be  the  reciprocal  of  the  thickness  of  the  ion  atmosphere 

where  this  thickness  is  the  distance  from  r = a^  to  a sphere  such  that  all  the 

other  ions  act  upon  J as  if  they  were  uniformly  distributed  upon  its  surface. 

Eqn.  (1-13)  can  be  split  into  the  contributions  to  0 , the  potential  at  a 

( 3 ) 

distance  r from  J,  of  the  central  ion  J and  of  the  ion  atmosphere  and,  at  r = 
( / 

a^,  the  potential  due  to  the  ion  atmosphere  is  found  to  be 


1 + Ka, 


(1-14) 


That  part  of  the  electrical  work  done  to  charge  J which  results  from  the  presence 
of  other  ions  is  then,  if  equated  to  the  electrical  free  energy, 


(el)  D 


1 Kqdq 


1 + Ka 


(1-15) 


where  qdq  represents  the  differential  charge  increment.  One  then  obtains, 
keeping  in  mind  that  K is  also  a function  of  the  charge  (l  owler  and  Guggenheim, 
1939,  p.  393;  Koryta  cl  al_.  , 1973,  p.  22) 


where 


2 

TtKa^)  = ■—  [ln{  1 + x)  - x + — ] with  x = Ka^  (1-17) 

x 

Note  that  A F , can  be  calculated  for  any  one  ion  in  solution  so  that  we  may 
(el) 

replace  the  subscript  J which  characterizes  a specific  central  ion  by  the  symbol 
i.  Furthermore,  the  free  energy  F rather  than  the  work  function  is  used 
because  the  small  electrostrictive  term  is  neglected  so  that  there  is  no  PdV 
work  considered. 

For  a component  i in  solution  one  has 

u = u ° + RT  In  f + RT  In  m.  (1-18) 

ri  l l i 

where  u and  u ° are  the  actual  and  the  standard  chemical  potentials,  f.  is  the 
i t i 

molal  activity  coefficient,  and  rm  is  the  molality.  Equivalent  expressions 
could  have  been  written  in  terms  of  the  rational  and  of  the  molar  activity  coef- 
ficients. 

If  the  departure  from  ideality  is  attributed  to  coulombic  interactions 


between  ions  then 


] -7 


dAF 


An 


LnUi 


(cl)i 


dn. 


11 T In  f. 


(1-19) 


dn.  represents  the  addition  of  an  infinitesimal  amount  of  i to  an  amount  of 
solution  which  contains  1 kg  of  water  and  could  have  been  represented  by  dm.. 
One  obtains  from  (1-16),  neglecting  the  electrostrictive  term, 


In  f 


(D-H) 


2 2 
z.  e 

i 

2D  k T 1 + Ka 
e D 


(1-20) 


which,  from  the  definition  of  It,  yields 


log  f 


(D-H) 


Az.V*5 

D i S 

T0.  5 

1 t BdVs 


(1-21) 


and 


log  f 


(D-H) 

± 


A z zllj'5 
D C Al  S 

, n .0.  5 

1 f DDaDIS 


d-22) 


where  f is  the  mean  activity  coefficient  of  an  electrolyte  while  z c and  z^  arc 
the  cationic  and  anionic  valences.  A ^ and  B ^ are  the  Debye  constants. 


In  molal  units 


(2.  3026  D k T) 
e B 


8tt  NAe  p 

BD  " ll 000  D k T 
' e B 


(1-24) 


AD  is  in  mole  kg  (deg  K)B  and  BD  in  cm  ^ mole  ^ kg^"  ^ (deg  K)^"  ^ . 

It 

The  Debye-Huckel  limiting  law  at  nearly  infinite  dilutions  can  be  obtained 

0.  5 

directly  from  Eqns.  ( 1 - 21 ) and  (1-22)  by  setting  B a I <<:  1 > or  can  l36 
derived  as  done  above  but  as  suming.  that  the  ions  are  far  enough  appart  so  that 
they  may  be  represented  by  point  charges.  The  law  then  has  the  forms 


, ,(D-H)  A 2.0.5 

!og  f.  = -Ad  z.  Is 


(1-25) 


. ( D-H)  . .0.5 

log  i±  = -Ad  zcza  Is 


(1-26) 


It  is  of  interest  to  students  of  this  field  to  have  available  to  them  the  values 
of  A^,  B^  and  of  f^B  ^ as  a function  of  for  various  molalities, which  we 
calculated  and  present  in  the  next  two  tables.  The  values  of  ^ will  also  be 
useful  in  what  follows. 


b)  Parameters  and  activity  coefficients  from  the  Dcbye-IIuckcl  equation 


Values  for  the  Debyc-lluckel  parameters  arc  presented  in  Table  1-1 


calculated  from  (1-26)  arc  presented  in  Table  1-2  and  1-3 


as  functions  of  the  ionic 


a misnomer  as  it 


c)  Interpolated  values  of  experimental  mean  activity  coefficients  of  1-]  electro 


lytes  (Culberson  equation 


Experimental  results  for  mean  activity  coefficients  will  be  needed  for 


comparative  purposes  with  various  theoretical  equations.  As  such  results  are 


too  widely  spaced  as  a function  of  ionic  strength,  C.  H.  Culberson  (Pytkowicz 


Atlas  and  Culberson,  1965)  applied  an  empirical  interpolation  equation  of 


to  several  1-1  and  2-1  electrolytes  at  25°  C 


Values  of  the  parameters  are  presented  in  Table  1-3  and  of  f 

±(interp) 

s f (experimental  values)  in  Table  1-4.  It  can  be  seen  that  one  i 


use  values  of  f to  represent  actual  values  as  will  be  done  in  the  next 

i(interp) 

section.  Four  decimals  places  are  reported  to  smooth  the  data  although  only 


the  first  three  are  known  because  f 


was  obtained  by  fitting  Eqn.  (1-27) 


Table  1-1.  Molal  values  of  the  Debye-Huckcl  parameters  and  11  as  a 
function  of  temperalur e. 


'(•C) 

ad 

[mole  * kg^  ‘ (< 

0 

0. 4918 

5 

0. 4953 

10 

0. 4989 

15 

0.  5027 

20 

0. 5067 

25 

0. 5108 

30 

0. 5151 

35 

0. 5196 

40 

0.  5243 

r - 1 , -0.  5 , 0. 5 , , „.0. 5. 

[cm  mole  kg  (deg  K) 


0.  3248 
0.3256 
0.  3264 
0. 3271 
0.3279 
0. 3287 
0.  3294 
0.  3302 
0. 3310 


* These  values  are  calculated  from  Eqns.  (1-23)  and  (1-24)  using  the  density  of 
pure  water  and  pertain,  therefore,  to  dilute  solutions. 
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Table  1-2.  Values  of  \ the  molal  activity  coefficient  from  the  Debye- 


Huckcl  equation,  calculated  in  this  work  for  1-1  electrolytes  at 


25°  C,  as  a function  of  the  ionic  strength  I and  of  a . 

j u 


a^  (distance  of  closest  approach) 


!s 

3.  80 

3.  82 

3.  84 

3.  86 

3.  88 

3.  90 

3.  92 

. 0001 

. 9885 

. 9885 

. 9885 

. 9885 

. 9885 

. 9885 

. 9885 

. 0005 

. 9747 

. 9747 

. 9747 

. 9747 

. 9748 

. 9748 

. 9748 

. 0010 

. 9649 

. 9649 

. 9649 

. 9649 

. 9649 

. 9649 

. 9649 

. 0015 

. 9575 

. 9575 

. 9575 

. 9575 

. 9575 

. 9575 

. 9575 

. 0020 

. 9514 

. 9514 

. 9514 

. 9514 

.9515 

. 9515 

.9515 

. 0025 

. 9462 

. 9462 

. 946  2 

. 9462 

. 946  2 

. 9462 

. 9462 

. 0030 

.9415 

.9415 

.9415 

.9415 

. 9416 

.9416 

. 9416 

. 0035 

. 9373 

. 9373 

. 9373 

. 9373 

. 9373 

. 9374 

. 9374 

. 0040 

. 9334 

. 9334 

. 9334 

. 9335 

. 9335 

. 9335 

. 9335 

.0045 

. 9298 

. 9298 

. 9298 

. 9299 

. 9299 

. 9299 

.9300 

. 0050 

.9264 

. 9265 

. 9265 

. 9265 

. 9265 

. 9266 

. 9266 

. 0055 

.9233 

. 9233 

. 9233 

. 9234 

. 9234 

. 9234 

. 9235 

. 0060 

. 9203 

. 9203 

. 9204 

. 9204 

. 9204 

. 9205 

. 9205 

. 0065 

.9175 

.9175 

.9175 

.9176 

.9176 

.9176 

.9177 

. 0070 

. 9148 

. 9148 

.9148 

.9149 

.9149 

.9150 

.9150 

. 0075 

.9122 

.9122 

. 91  23 

.9123 

.9124 

.9124 

.9124 

. 0080 

. 9097 

. 9098 

. 9098 

. 9098 

. 9099 

. 9099 

.9100 

. 0085 

. 9073 

. 9074 

. 9074 

. 9075 

. 9075 

. 9076 

. 9076 

. 0090 

. 905  1 

. 9051 

. 9052 

. 9052 

. 9053 

. 9053 

. 9054 

. 0095 

. 9029 

. 9029 

. 9030 

. 9030 

. 903  1 

. 9031 

. 9032 

.0100 

. 9007 

. 9008 

. 9008 

. 9009 

. 9009 

. 9010 

.9011 

.0110 

. 8967 

. 8967 

. 8968 

. 8968 

. 8969 

. 8970 

. 8970 

.0120 

. 8929 

. 8929 

. 8930 

. 8930 

. 893  1 

. 8932 

. 8932 

.0130 

. 8892 

. 8893 

. 8894 

. 8894 

. 8895 

. 8896 

. 8897 

.0140 

. 8858 

. 8859 

. 8860 

. 8860 

. 8861 

. 8862 

. 8863 

.0150 

. 8826 

. 8826 

. 8827 

. 8828 

. 8829 

. 8829 

. 8830 

.0160 

. 8794 

. 8795 

. 8796 

. 8797 

. 8798 

. 8798 

. 8799 

.0170 

. 8764 

. 8765 

. 8766 

. 8767 

. 8768 

. 8769 

. 8770 

.0180 

. 8736 

. 8737 

. 8738 

. 8739 

. 8739 

. 8740 

. 874  1 

.0190 

. 8708 

. 8709 

.8710 

.8711 

.8712 

.8713 

. 8714 

. 0200 

. 868  2 

. 8683 

. 8684 

. 8685 

. 8686 

. 8687 

. 8688 

. 0300 

. 8458 

. 8459 

. 84  6 1 

. 8462 

. 84  63 

. 8465 

. 8466 

. 0400 

. 8284 

. 8286 

. 8288 

. 8289 

. 8291 

. 8293 

. 8294 

. 0500 

. 8142 

. 8144 

. 8146 

.8147 

.8149 

.8151 

.8153 
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. 0600 

. 8020 

. 8023 

. 8025 

. 80  27 

. 8029 

. 8031 

. 8033 

.07  00 

.7914 

.7917 

.7919 

. 7922 

. 7924 

. 7927 

. 7929 

. 0800 

. 7821 

.7823 

. 7826 

. 7828 

. 783  1 

. 7834 

. 7836 

.0900 

. 7736 

. 7739 

. 774  2 

. 7745 

. 7748 

. 7750 

. 7753 

. 1000 

. 7660 

. 7663 

. 7666 

. 7669 

. 7672 

. 7675 

. 7678 

. 2000 

.7136 

. 7140 

.7145 

.7149 

.7154 

.7)58 

.7163 

. 3000 

. 6821 

.6827 

. 6833 

.6838 

.6844 

.6849 

.6855 

. 4000 

. 6600 

.6606 

.6612 

.6619 

.6625 

. 663  1 

. 6637 

. 5000 

.6430 

. 6437 

. 6444 

.6451 

. 6458 

. 646  5 

.6471 

.6000 

.6294 

. 6301 

. 6309 

.6316 

.6324 

. 633  1 

. 6338 

. 7000 

.6180 

.6188 

.6196 

.6204 

.6212 

.6  220 

.6228 

. 8000 

.6084 

.6093 

.6101 

.6109 

.6118 

.6126 

.6134 

. 9000 

.6001 

.6010 

.6018 

.6027 

. 6036 

.6044 

.6053 

1.0000 

. 5928 

. 5937 

. 5946 

. 5955 
aD 

. 5964 

. 5972 

. 5981 

xs 

3.  94 

3.  96 

3.  98 

4.  00 

4.  20 

4.40 

4.60 

. 0001 

. 9885 

. 9885 

. 9885 

. 9885 

.9885 

. 9885 

. 9885 

. 0005 

. 9748 

. 9748 

. 9748 

.9748 

. 9748 

. 9748 

.9749 

. 0010 

. 9649 

. 9649 

. 9649 

. 9649 

. 9650 

.965  1 

. 9651 

.0015 

.9576 

. 9576 

. 9576 

. 9576 

. 9577 

. 9578 

. 9579 

. 0020 

. 9515 

. 9515 

.9515 

. 9515 

. 9517 

.9518 

.9519 

. 0025 

. 9463 

. 9463 

. 9463 

. 9463 

. 9465 

. 9466 

. 9468 

. 0030 

. 9416 

.9416 

. 9417 

. 9417 

.9419 

. 9421 

. 9422 

. 0035 

. 9374 

. 9374 

. 9375 

. 9375 

. 9377 

. 9379 

. 9381 

. 0040 

. 9336 

. 9336 

. 9336 

. 9336 

. 9339 

. 9341 

. 9344 

. 0045 

. 9300 

. 9300 

. 9300 

. 9301 

. 9303 

. 9306 

. 9309 

. 0050 

. 9266 

. 9267 

. 9267 

. 9267 

. 9270 

.9273 

. 9276 

, 0055 

. 9235 

. 9235 

. 9236 

. 9236 

.9239 

. 9242 

. 9246 

. 0060 

. 9205 

. 9206 

.9  206 

. 9 206 

. 9210 

. 9213 

. 9217 

.0065 

.9177 

.9178 

.9178 

.9178 

.9182 

.9186 

.9190 

. 0070 

.9150 

.9151 

.9151 

.9152 

.9156 

.9160 

.9164 

. 0075 

.9125 

.9125 

.9126 

.9126 

.9130 

.9135 

.9139 

. 0080 

.9100 

.9101 

.9101 

.9102 

.9106 

. 91  1 i 

.9115 

. 0085 

. 9077 

. 9077 

. 9078 

. 9078 

. 9083 

. 9088 

. 9092 

. 0090 

. 9054 

. 9055 

. 9055 

. 9056 

. 906  1 

. 9065 

. 9070 

. 0095 

. 903  2 

. 9033 

. 9033 

. 9034 

. 9039 

. 9044 

. 9049 

.0100 

.9011 

.9012 

.9012 

.9013 

. 9018 

. 90  23 

. 9029 

.0110 

. 8971 

. 8971 

. 897  2 

. 8973 

. 8978 

. 8984 

. 8990 

.0120 

. 8933 

. 8934 

. 8934 

. 8935 

. 894  1 

. 8947 

. 8954 

.0130 

. 8897 

. 8898 

. 8899 

. 8899 

. 8906 

.8913 

.8919 

.0140 

. 8863 

. 8864 

. 8865 

. 8865 

. 8873 

. 8880 

. 8887 

.0150 

. 8831 

. 8832 

. 8832 

. 8833 

. 8841 

. 8848 

. 8856 

.0160 

. 8800 

. 8801 

. 8802 

. 8802 

. 8810 

.8818 

. 8826 

.0170 

. 8770 

. 8771 

. 8772 

. 8773 

. 8781 

. 8790 

. 8798 

.0180 

. 874  2 

. 8743 

. 8744 

. 8745 

. 8753 

. 876  2 

. 8771 

.0190 

.8715 

. 8716 

.8717 

.8718 

. 8727 

. 8736 

. 8744 

. 0200 

. 86  88 

. 8689 

. 8690 

. 8691 

. 8701 

.8710 

.8719 

. 0300 

. 8467 

. 8468 

. 8470 

.8471 

. 8484 

. 84  97 

. 8509 

. 0400 

. 8296 

. 8297 

. 8299 

. 8301 

.83  17 

. 8332 

. 8348 

. 0500 

. 8155 

.8157 

. 8159 

.8161 

.8179 

. 8198 

. 8216 

. 0600 

. 8035 

. 8038 

. 8040 

. 804  2 

. 8063 

. 8084 

. 8104 

. 0700 

. 793  1 

. 7934 

. 7936 

. 7938 

. 796  2 

. 7985 

. 8007 

. 0800 

. 7839 

. 7842 

. 7844 

. 7847 

. 7872 

. 7897 

. 7921 

. 0900 

. 7756 

. 7759 

. 7762 

. 7764 

. 7792 

. 7819 

. 7845 

. 1000 

. 7681 

. 7684 

. 7687 

. 7690 

. 7719 

. 7748 

. 7775 

. 2000 

. 7 167 

.7172 

.7176 

.7180 

.7224 

. 7266 

. 7307 

. 3000 

. 6860 

. 6865 

. 6871 

.6876 

. 6929 

. 6981 

. 7030 

. 4000 

. 6644 

. 6650 

. 6656 

.6662 

. 6722 

.6781 

.6837 

. 5000 

.6478 

. 6485 

. 6492 

.6499 

. 6565 

. 6629 

. 6690 

. 6000 

.6346 

. 6353 

. 6360 

. 6368 

. 6439 

. 6507 

. 6573 

. 7000 

. 6236 

.6243 

. 6251 

.6259 

. 6334 

. 6407 

. 6476 

. 8000 

.6142 

.6150 

.6158 

.6166 

. 6245 

. 6321 

. 6394 

. 9000 

.6061 

.6070 

.6078 

.6087 

.6169 

. 6248 

. 6323 

1. 0000 

. 5990 

. 5999 

. 6008 

. 6016 
aD 

.6101 

.6183 

.6261 

!s 

4.  80 

5.  00 

5.  20 

5. 40 

5. 60 

5.  80 

6.  00 

. 0001 

. 9885 

. 9885 

. 9885 

. 9885 

. 9885 

. 9885 

. 9885 

.0005 

. 9749 

. 9750 

. 9750 

. 9750 

. 9751 

. 9751 

. 975  1 

. 0010 

. 9652 

. 9653 

. 9653 

. 9654 

. 9655 

. 9655 

. 9656 

.0015 

. 9580 

. 9581 

. 9582 

. 9583 

. 9584 

. 9585 

. 9586 

. 0020 

. 9521 

. 9522 

. 9523 

. 95  24 

. 9526 

. 9527 

. 9528 

. 0025 

. 9470 

. 9471 

. 9473 

. 9474 

. 9476 

. 9477 

. 9479 

.0030 

. 9424 

. 9426 

. 9428 

. 9430 

. 9432 

. 9433 

. 9435 

. 0035 

. 9383 

. 9386 

. 9388 

. 9390 

. 9392 

. 9394 

. 9396 

. 0040 

. 9346 

. 9348 

. 9351 

. 9353 

. 9355 

. 9358 

. 9360 

. 0045 

.9311 

.9314 

. 9317 

.9319 

. 9322 

. 9324 

. 9327 

. 0050 

. 9279 

. 9282 

. 9285 

. 9288 

. 9291 

. 9293 

. 9296 

. 0055 

.9249 

. 9252 

. 9255 

. 9258 

. 926  1 

. 9264 

. 9267 

. 0060 

. 9220 

. 9224 

. 9227 

. 9230 

. 9234 

. 9237 

. 9240 

O 

o 

LTl 

.9193 

.9197 

. 9200 

. 9204 

. 9207 

. 921  1 

.9214 

. 0070 

.9167 

.9171 

.9175 

.9179 

.9183 

.9186 

.9190 

. 0075 

.9143 

.9147 

.9151 

.9155 

.9159 

.9163 

.9167 

. 0080 

.9119 

.9124 

.9128 

.9132 

.9136 

.9140 

.9145 

. 0085 

. 9097 

. 9 J 0 J 

. 9106 

.9110 

.9115 

.9119 

.9123 

. 0090 

. 9075 

. 9080 

. 9085 

. 9089 

. 9094 

. 9098 

.9103 

. 0095 

. 9054 

.9059 

. 9064 

. 9069 

. 9074 

. 9079 

. 9083 
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0100 

.9034 

. 9039 

.9044 

. 9049 

. 9054 

. 9059 

. 9064 

01  10 

. 8996 

. 9001 

. 9007 

.9012 

.9018 

. 90  23 

. 9028 

0120 

. 8960 

. 8966 

. 8972 

. 8977 

. 8983 

. 8989 

. 8995 

0130 

. 89  Z6 

. 8932 

. 8938 

. 8945 

. 8951 

. 8957 

. 8963 

0140 

. 8893 

. 8900 

. 8907 

. 8914 

. 8920 

. 8927 

. 893  3 

0150 

. 8863 

. 8870 

. 8877 

. 8884 

. 8891 

. 8898 

. 8905 

0160 

. 8834 

. 884  1 

. 884  9 

. 8856 

. 8863 

. 8870 

. 8877 

0170 

. 8806 

. 8814 

. 8821 

. 8829 

. 8837 

. 8844 

. 8852 

0180 

. 8779 

. -7  87 

. 8795 

. 8803 

.8811 

. 8819 

. 8827 

0190 

. 8753 

. 876  2 

. 8770 

. 8779 

. 8787 

. 8795 

. 8803 

0200 

. 8728 

. 8737 

. 8746 

. 8755 

. 8764 

. 8772 

. 8780 

0300 

. 8521 

. 8534 

. 8545 

. 8557 

. 8569 

. 8580 

. 8591 

0400 

. 8363 

. 8377 

. 8392 

. 8406 

. 8420 

. 8434 

. 8448 

0500 

. 8233 

. 8250 

. 8267 

. 8284 

. 8300 

. 8316 

. 8332 

0600 

.8124 

.8143 

.8162 

. 8181 

. 8199 

. 8217 

. 8234 

0700 

. 8029 

. 8050 

. 8071 

. 8092 

.8112 

. 8131 

.8151 

0800 

. 7945 

. 7968 

.7991 

.8013 

. 8035 

. 8056 

. 8077 

0900 

. 7870 

. 7895 

. 7920 

. 7943 

. 7967 

. 7989 

.8012 

1000 

. 7803 

. 7829 

. 7855 

. 7880 

.7905 

.7929 

. 7953 

2000 

. 7346 

. 7385 

. 7422 

. 7458 

. 7494 

. 7528 

. 7562 

3000 

. 7078 

.7125 

.7170 

. 7213 

.7256 

. 7297 

. 7337 

4000 

.6891 

. 6944 

.6995 

. 7044 

. 7091 

.7137 

.7182 

5000 

.6750 

.6807 

.6862 

.6916 

.6967 

.7017 

. 7066 

6000 

.6637 

.6698 

. 6757 

.6814 

.6869 

.6922 

.6974 

7000 

. 6543 

.6608 

.6670 

. 6730 

.6788 

.6844 

.6898 

8000 

.6464 

.6532 

. 6597 

.6659 

.6720 

.6778 

.6834 

9000 

.6396 

. 6466 

. 6534 

.6598 

. 666 1 

.6721 

.6780 

0000 

. 6336 

.6409 

.6478 

.6545 

.6610 

.6672 

.6732 
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Tabic  1-3.  Coefficients 

for  Eqn.  (1-27)  and  the 

standard  deviation  of  fit 

relative  to  measured  mean  activity  coefficients  at  25°  C. 

KC1 

NaCl 

Ac 

-0. 5108 

-0. 5108 

Bc 

1. 3070 

1. 3500 

o 

o 

0 

0. 043663 

Dc 

0 

-0. 0093698 

Ec 

0. 0020748 

0 

Molality  range 

0-1.0 

0 

1 

o 

Std.  deviation  of  fit 

0. 08%  in  f 

0.  09%  in  f 
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Tabic  1 -4. 


Values  of  f obtained  from  Eqn.  (1-27)  and  experimental 

±(  inter  p)  ~ ' 

values  f from  the  literature.  Data  on  the  rnolal  scale  at  25°  C. 

±(cxp) 

KC1  NaCl 

f . f f f 

±(inter  p) i.(  exp) jointer  p) l(cxp) 


0. 0001 

0. 9885 

0.0005 

0.  9748 

0.  975 

0.  9748 

0.  975 

0.  001 

0. 9649 

0.  965 

0. 9651 

0.  965 

0.  002 

0. 95  15 

0.  952 

0. 9518 

0.  952 

0.  003 

0. 94  17 

0.  9421 

0.  004 

0. 9336 

0. 9341 

0.005 

0.  9267 

0.  928 

0. 9273 

0.  927 

0.  006 

0.  9206 

0.  9214 

0.  007 

0. 9151 

0. 9160 

0.  008 

0.9101 

0. 91  1 1 

0.  00? 

0. 9055 

0. 9066 

0.  010 

0. 9012 

0.  902 

0.  9024 

0.  902 

0.  020 

0. 8690 

0. 8713 

0.  030 

0.  8470 

0. 8503 

0.  040 

0. 8299 

0. 8341 

0.  050 

0. 8159 

0. 8210 

0.  060 

0. 8040 

0. 8099 

0.  070 

0. 7936 

0. 8004 

0.  080 

0. 7844 

0.  7920 

0.  090 

0. 7761 

0.  7845 

0.  10 

0.  7686 

0.  770 

0. 7778 

0.  778 

0.  15 

0. 7391 

0. 75  18 

0.  20 

0.  7177 

0.718 

0.  7336 

0.  735 

0.  25 

0. 7009 

0.  7199 

0.  30 

0.  6873 

0.688 

0.  7091 

0.  710 

0.  35 

0.  6758 

0.  7004 

0.40 

0. 6660 

0.  666 

0.  6931 

0.693 

0.45 

0. 6574 

0.  6871 

0.  50 

0. 6498 

0.  649 

0.6819 

0.681 

0.  55 

0. 6431 

0.  6775 

0.  60 

0.  6370 

0.  637 

0.6736 

0.  673 

0.  65 

0.6315 

0. 6703 

0.  70 

0.  6264 

0.626 

0. 6674 

0.  667 

0.  75 

0.  6218 

0. 6649 

0.  80 

0.  6176 

0.618 

0 . 6 6 26 

0.  662 

0.  85 

0.6137 

0. 6606 

0.  90 

0.6100 

0.  610 

0.6589 

0.  659 

0.  95 

0. 6066 

0. 6574 

1.0 

0. 6035 

0.  604 

0. 6560 

0.657 

2.  0 

0.  573 

0.  668 

3.  0 

0.  569 

0.  714 

4.0 

0.  577 

0.  783 

5.  0 

0.  874 

6.  0 

0.  986 

— 


... 


to  experimental  data  reported  to  three  decimals 


d)  Comparison  of  Debye- 1 luckel  and  experimental  results 


It  is  of  interest  to  compare  results  obtained  from  Ecjn.  (1-26)  and  shown 


Guntelberg  (1926)  suggested  a simplified  form  of  the  Debye-Huckel 


equation,  namely 


3.04  A at  25 ° C to  be  applied  to  all  electrolytes 


which  corresponds  to  a 


This  form  was  adopted  by  Guggenheim  ( 1935)  and  by  Davies  (1938)  among 


others  as  the  first  term,  which  describes  non-specific  coulombic  interactions 


in  their  equations.  Values  of  log  f 


for  1-1  electrolytes  are  presented  in 


Table  1-6.  If  one  remembers  that  a 


ion  with  its  counter-ion,  a virtue  of  Eqn.  (1-28)  is  that,  for  empirical  fitting 


of  data  by  means  of  additional  specific  interaction  terms  in  mixed  electrolyte 


solutions,  one  docs  avoid  the  uncertain  meaning  of  a 


This  quantity  is  only 


well  defined  in  the  Debyc-IIuckel  equation  for  a single  electrolyte  at  a time 


The  results  in  Table  1-5  indicate  that  the  Debye-Huckel  equation  fits 


KC1  data  quite  well;  to  within  less  than  0.  1%  at  m = 0.  1 and  to  within  less  than 


1%  at  m = 1.0.  The  agreement  is  not  as  good  for  NaCl  and  becomes  worse 
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II 

Table  1-5. 

A comparison  of  values 

calculated  from  the  Debye-Huckel 

equation 

with  exper 

imental  ones 

as  represented  by  Eqn. 

(1-27) 

repre- 

sents  the  percent  deviat 

, .(El- 
ion of  f 

± 

H)  r 

f rom  f , . 

±(  inter  p) 

KCl  (a 

o 

= 3. 94  A) 

NaCl  (a 

o 

= 4. 60  A) 

m 

^£(  inter p) 

f(D-II) 

A% 

^i(interp) 

{( D-H) 
± 

L°k 

0.  001 

0.  9649 

0. 9649 

0.  00 

0. 9651 

0. 9651 

0.  00 

0.  005 

0.  9267 

0.  9266 

0.  01 

0. 9273 

0. 9276 

0.  03 

0.  01 

0. 9012 

0. 901 1 

0.  01 

0.  9024 

0. 9029 

0.  06 

0.  05 

0.  9159 

0. 8155 

0.  05 

0. 8210 

0. 8216 

0.  07 

0.  1 

0.  7686 

0.  7681 

0.  07 

0. 7778 

0. 7775 

0.  04 

0.  2 

0. 7177 

0. 7167 

0.  14 

0. 7336 

0. 7307 

0.41 

0.  3 

0.  6873 

0.  6860 

0.  18 

0. 7091 

0. 7030 

0.  86 

0.4 

0. 6660 

0.  6644 

0.  24 

0.6931 

0.  6837 

1.  36 

0.  5 

0.6498 

0.  6478 

0.  31 

0.  6819 

0. 6690 

1.  89 

1.0 

0.6035 

0.  5990 

0.  75 

0. 6560 

0.6216 

5.  24 

(G)  " 

Table  1-6.  Values  of  f obtained  from  the  Guntclberg  equation  for  1-1 


electrolytes  at  25°  C. 


0.  001 


0.  005 


0. 9646 


0.  9253 


0. 8986 


0. 8066 


0. 7538 
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1 


foi*  other  electrolytes,  especially  multivalent  ones.  1 he  Guntclborg  ctju^tion, 
as  is  to  he  expected  provides  a poorer  fit. 

It  is  our  contention  that,  as  will  be  discussed  in  the  next  section,  the 
Debye-IIuckel  equation  is  not  theoretically  satisfactory  and  that  it,  as  well  as 
extended  forms  based  upon  it,  provide  activity  coefficients  that  agree  with 
measured  ones  through  empirical  curve  fitting. 

II 

e)  Critique  of  the  Dcbye-Huckcl  theory 

II 

Several  conceptual  problems  of  the  Debye-Ruckel  theory  arise  from  the 

use  of  the  Boltzman  distribution  for  the  distribution  of  ions  around  a central 

ion.  The  Boltzman  law  is  derived  for  non-interacting  particles  or  for  canonical 

ensembles  of  systems  which  may  contain  interacting  particles.  And  yet  the  law 

is  applied  to  individual  particles  in  a problem  the  very  essence  of  which  is,  for 

intermediate  concentrations,  to  calculate  activity  coefficients  resulting  from 

sizeable  coulombic  interactions.  As  an  example,  at  m = 0.  1 the  nearest 

14 

neighbor  energy  for  a 1-1  electrolyte  is  1.459  X 10  ergs,  a significant 

14 

fraction  of  k T = 5.  836  X 10  ergs. 

13 

Some  problems  that  arise  from  the  Boltzman  law  are  resolved  by  the 
linearization  of  the  density.  This  is  especially  true  when  the  linearization  is 
used  to  render  the  density  of  charge  and  the  resulting  potential  compatible  with 
the  principle  of  superposition  of  charges.  This  procedure  depends,  however, 


upon  the  condition 


1-2) 


especially  in  the  case  of  asymmetric  electrolytes  for  which  terms  in  even 
powers  do  not  vanish  in  principle  in  Eqn.  (1-7).  E is  the  electrical  potential 
energy  of  an  i-ion. 

It  is  easy  to  show  that  condition  (1-29)  is  not  met  even  in  what  is  usually 

II 

considered  the  range  of  validity  of  the  extended  Debye-Huckel  equation.  Con- 
sider the  energy  of  interaction  of  a single  cation  with  an  anion 


E = 


2 

c 


D r 
e o 


(1-30) 


for  a 1-1  electrolyte,  r^  is  the  nearest  neighbor  distance.  This  energy  is 

smaller  than  that  of  an  ion  in  an  actual  solution  as  will  become  apparent  when 

we  develop  the  partial  order  theory  in  later  chapters.  We  do  not  use  here 

Eqn.  (1-2°)  because  0 given  by  Eqn.  (1-13)  already  presupposes  the  validity 

of  condition  (1-29).  Values  of  E are  compared  to  k T in  Table  1-7  and  it  can 

B 

be  seen  that  the  neglect  of  the  second  power  term  is  not  justified  for  asymmetric 
electrolytes  even  at  m = 0.01  and  the  neglect  of  the  third  power  term  for  sym- 
metric  electrolytes  at  m > 0.  05.  The  errors  would  of  course  be  more  severe 
for  higher  valence  salts. 

ii 

An  important  contradiction  of  the  Debye-Iluckel  theory  which  is  applicable 
even  to  the  limiting  law  arises  in  the  ion  atmosphere  concept.  Let  us  examine 
Eqn.  (1-2)  or,  rather,  its  continuous  linearized  counterpart.  After  linearization 
it  can  be  shown  that 
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Table  1-7.  Values  of  powers  of  x = E/k^T  and  of  their  percent  ratios  to 
E/k  T versus  the  molality  for  1-1  electrolytes  at  25°  C. 


m 

r x 108 

O 

(cm) 

E x 1014 
(ergs) 

E/kBT 

2 

X 

3 

X 

2. 

X /x 
(%) 

x 3 /x 
(%) 

0.  01 

43.  73 

0.  678 

0.  165 

0.  027 

0. 0045 

16.  5 

2.  7 

0.  05 

25.  58 

1.  159 

0.  282 

0.  079 

0.  022 

28.  2 

7.  9 

0.  1 

20.  31 

1. 459 

0.  354 

0.  125 

0.  044 

35.  4 

J 2.  5 

0.  2 

16.  13 

1.  837 

C.  446 

0.  199 

0.  089 

44.6 

19.  9 

0.  5 

11.91 

2.  489 

0.  605 

0.366 

0.  221 

60.  5 

36.  5 
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D K 

P(J)  = 0(J) 


d-31) 


and,  introducing  Eqn.  (1-13), 


Z,e  2 *aD  -*r(J| 
J K c e 

’(J)  = _4IT  "]  +KaD  r 


(1-32) 


One  observes  that  p , the  net  charge  density,  decays  monotonically  from  a 

( J ) 

maximum  value  at  r = a to  zero  at  r = <*>.  This  is  simply  a formal  state- 

( J ) U l J / 

ment  of  the  ion  atmosphere  concept,  namely,  that  the  surroundings  of  J contain 
an  excess  of  counter-ions  and  that  this  excess  decays  with  distance.  Thi s 
concept  leads  to  a serious  condradiction  because  if  now  we  select  one  of  the 
counter-ions  near  the  original  J as  the  new  central  ion  J,  it  is  surrounded  by 

it 

an  excess  of  ions  like  it  in  contradiction  to  the  basic  tenet  of  the  Debye-Huckel 
theory.  Thus,  there  is  undue  bunching  of  like  ions. 

1 1 

Another  conceptual  difficulty  is  that  the  Debye-Huckel  theory  leads  to 
potentials  which  are  not  uniquely  defined  as  a function  of  position  in  a solution. 
In  effect,  let  us  retain  a volume  element  dV,  at  which  we  calculate  by  means  of 
Eqn.  (1-13)  the  potential  0 due  to  all  the  ions  in  solution,  but  let  us  take  two 

( J ) 

ions  of  the  same  kind  as  successively  the  central  ion.  Let  these  two  ions  be 
J1  and  J2,  placed  at  different  distances  from  dV  (Fig.  1-1).  The  potentials  at 


J due  to  the  solution  will  be 


r 


*a 

z T e D 
J e 


J1_De1+KaD  r(Jl) 


(1-33) 


J2  D 1 + Ka 
e D 


(1-34) 


that  the  potential  at  dV  due  to  all  the  ions  present  in  the  solution  depends  upon 


the  choice  of  the  reference  ion. 


These  problems  arise  because  Debye  and  Huckel  did  not  take  into  con- 
sideration the  quasi-lattice  nature  of  ionic  solutions,  that  is,  that  there  is  some 
degree  of  regularity  in  the  spacial  distribution  of  ions  which  precludes  excessive 
bunching  of  counter-ions  in  the  ionic  atmosphere  as  well  as  a decay  of  the  net 
charge  density  towards  zero  at  large  distances  from  any  given  ion. 


If  we  look  at  the  essence  of  the  Debye-Huckel  theory,  it  consists  in 


obtaining  a quantity 


V 1 


D r. 
e la 


(1-35) 


where  r.  , the  radius  of  the  ion  atmosphere,  is  a + 1/K  and  integrating  it  to 


obtain  AF  . If  wc  consider  Eqn.  (1-35)  alone,  where  r is  the  radius  of  a 
(el) 


spher e 


such  that  -z  £ distributed  uniformly  over  it  simulates  the  effect  of  the 
J 


1 
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solution  potential  upon  .7,  without  concern  about  how  r.  was  obtained,  then 

Eqn.  (1-35)  is  logically  sound.  The  difficulty  arises  with  a proper  evaluation 

of  r which  takes  into  account  the  actual  and  partially  regular  distribution  of 
ia 

II 

cations  and  anions  in  the  solution.  The  apparent  success  of  the  Debye-Huckel 

equation  indicates  that  their  choice  of  r.  = a_  + 1/K  yields,  by  means  of  the 

la  D 

adjustable  parameter  a^,  values  of  5AF,  ,,/dm  which  are  correct  at  low  concen- 
J 1 D (el) 

trations. 

/ 

Note  that  the  above  conclusion  does  not  imply  that  the  values  of  0 and 

aD 

of  AF  are  correct  but  only  that  the  derivative  of  the  latter  yields  proper 
(el) 

values,  that  is,  that  the  slope  of  AF  versus  m is  the  proper  one.  We  will 

\ ) 

7 

obtain  from  the  partial  order  theory  at  m = 0.  01  a value  AF  = -3.  508  X 10 

(el) 

" 7 

ergs  while  from  the  Debye-Huckel  equations  one  obtains  AF  ^ = -4.537  x 10 

c 

ergs  and  yet  both  approaches  yield  f = 0.  902  (if  a^  = 4.  60  A is  used  in  the 

ii 

Debye-Huckel  equation).  These  results  confirm  the  opening  statement  of  this 

paragraph,  if  our  value  of  AF  is  accepted  as  the  correct  one,  a point  which 

\ ) 

the  reader  will  be  able  to  evaluate  later  on  in  this  work. 

ri 

In  the  limiting  case  the  Debye-Huckel  theory  should  be  valid  because  the 
Boltzman  distribution  should  be  nearly  correct  for  weak  interactions  and  our 
objections,  which  result  from  the  use  of  this  distribution,  should  correspond 

fl 

to  slight  erro-s.  Furthermore,  the  static  nature  of  the  Debye-Huckel  theory 
is  acceptable  in  the  limiting  case  because  the  effect  of  fluctuations,  as  seen 


by  a central  ion,  are  unimportant  at  high  dilutions  (Frank  and  Thompson,  1959). 


2-  J 


2-  Original  Latin  r Theory;  Perfect  Long-Range  Order  (PLRO) 

In  the  original  lattic  e theory  (e.  g.  , Robinson  and  Stokes,  1 9S9 , p.  226) 
the  ions  arc  assumed  to  form  a regular  lattice,  with  perfect  long-range  order, 
in  the  solution.  This  theory  does  not  work,  that  is,  it  does  not  yield  activity 
coefficients  compatible  with  experimental  ones.  Still,  it  will  be  described 
because  it  serves  as  the  basis  for  our  partial  long-range  order  (pLRO)  which 
does  work. 

Notice  that  the  '’LRO  theory  contains  two  assumptions;  the  regular  dis- 
tribution of  ions  in  equidistant  sites  and  the  perfect  long-range  order  and, 

furthermore,  that  the  first  assumption  is  compatible  with  an  intrinsic  one 
" " 

made  by  Debye-Huckel.  Indeed,  Debye-Huckel  assumed  stationary  charges, 
which  must  therefore  occupy  average  regularly  spaced  positions,  because  the 
stationary  condition  is  a necessary  one  when  the  Poisson  equation  is  used. 

The  PLRO  theory  is  based  upon  that  of  ionic  crystal  lattices  which  will 
be  briefly  recapitulated  next, 
a)  Energy  of  interaction  in  an  ionic  crystal 

For  an  ionic  crystal  in  which  cations  and  anions  have  the  same  valences, 
the  energy  of  interaction  for  the  whole  lattice  is  (e.g.  Kittel,  1959;  Pauling, 
1940,  p.  337) 
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(2-1) 


2-2 


where  N is  the  number  of  ions  of  each  kind  in  the  volume  V under  consideration, 

A is  the  Madelung  constant,  and  B is  the  Born  constant.  The  last  term  is 
M B 

the  Born  repulsion  potential. 

Madelung' s constant  (Madelung,  1918)  represents  the  effects  of  interactions 

of  nearest  neighbors,  next  nearest  neighbors,  and  so  on  upon  the  lattice  energy. 

The  general  shapes  of  E . is  sketched  in  Fig.  2-1.  A*  r = r binding 

b r (xtal)  ° 

energy  passes  through  a minimum  and  yields  a stable  lattice, 
b ) Ionic  solution  with  perfect  long-range  order 

Let  us  divide  a solution  into  equal  cubes,  place  an  ion  at  the  center  of 
each  cube,  and  alternate  cations  and  anions.  The  resulting  ionic  lattice  is 
represented  by  Fig.  2-2.  The  structure  corresponds  to  a face- center ed  cubic 
lattice  with  anions  distributed  hexagonally  around  cations  and  vice-versa.  As 
the  short-range  repulsive  forces  would  only  be  significant  at  extremely  high 
concentrations,  the  counterpart  of  Eqn.  (2-1)  is  now 


m 2 2a 

N z e A.. 

„ A M 

AE  = c 

l Dr 

e o 


(2-2) 


where  AE  is  the  lattice  (configurational  coulombic)  energy  relative  to  that  at 
infinite  dilution  energy  and  c is  the  molarity.  Thus,  N^c  represents  the  number 
of  ior.s  of  each  kind  present  per  liter  of  solution.  Note  that  AF^  and  E^  can  be 
used  interchangeably  because  E^at  c = 0 is  zero. 

As  r is 
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Eqn.  (2-2)  becomes 


,1/3  4/3  2 2. 

,r  2 A Z G M 4/3  . 2 4/3 

hTd C = *BlV  C 


(2-4) 


with 


21/3N  4/V 

A 

10  D 

e 


(2-5) 


We  use  molalities  throughout  this  work  to  avoid  changes  in  concentration 
due  to  pressure  and  temperature  effects  so  let  us  convert  Eqn.  (2-4)  to  a molal 
base.  Then,  we  start  from 


M 2 2. 

N z e A 
A M 
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As 


pm 

1 + 0.  001  mM2 
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(2-7) 
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where  M .s  the  molecular  weight  of  the  solute  and  p is  the  density  of  the  solu- 

z 

tion,  r is  now' 
o 


r 

o 


1000  + mM 
2 

2N  mp 
A 


1/3 


(2-8) 


and 


or 


10B/2AMP,/3  4/3 

AE  - 1/3  m 

(1000  + mM  ) 


(2-9) 


P zZA  d1/3 

I M° 4/3 

AE  - - ] / 3 m 

( 1 + 0.  001  mM) 


(2-10) 


In  dilute  solutions  we  may  write  Eqn.  (2-8)  in  the  form 


AE 


» ..2a  . l/3_.4/3 


AM^o  ™ 


(2-11) 


and  calculate 


8AE, 


log  f = 


1 1 


£ 2.3026  HT  8m 


X 0.  5 


(2-12) 
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No  configurational  entropy  change  i.  used  in  this  model  because,  as  the  degree 
of  order  does  not  change  with  concentration,  Ah  e ' |,olc 


electrical  energy 
can  be  equated  to 
are  neglected. 


change  and  is,  therefore,  the  change  in  internal  energy.  It 
the  enthalpy  change,  because  electrostrictive  effects 


As  A = 1.7476  for  face- center ed  cubic  lattices,  one  obtains  from  Eqns 
M 

(2-11)  and  (2-12) 


log  f± 


.a  1/3 

2BlZ  AMPo  ml/3 

3 X 2.  3026  RT 


-0. 3842  m 


1/3 


(2-13) 


for  a 1-1  electrolyte.  In  molar  units 


!°g  f± 


-0. 3846 


1/3 

c 


(2-14) 


These  results  differ  from  those  reported  in  the  literature,  e.  g.  (Robinson 
and  Stokes,  1959,  p.  226) 


log 


-0. 29c 


1/3 


(2-15) 


because  the  latter  equation  is  improperly  derived. 


In  effect,  it  is  obtained  by 


setting 


by  stating  that  A E^,  represents  the  configurational  energy  when  one  mole  (N^ 
ions  of  each  kind)  of  salt  is  dissolved  in  an  infinite  amount  of  water  (i.  e.  , AE^ 
is  the  chemical  potential  A|^>).  And  yet,  r^  is  calculated  from  Eqn.  (2-3)  for 
a finite  molarity  c so  that  the  derivation  of  Eqn.  (2-15)  is  not  self-consistent. 

It  i s of  interest  to  compare  results  obtained  from  Eqns.  (2-14)  and  (2-15) 
with  experimental  data  noting  that,  to  a first  approximation,  we  do  not  have  to 
distinguish  c and  m for  rough  comparisons  at  high  dilutions.  This  is  done  in 
Table  2-1.  Neither  equation  agrees  well  with  experiments  as  is  to  be  expected 
because  the  energy  of  interaction  is  small  in  dilute  solutions  and  should  be  easily 
overcome  by  the  thermal  energy  k^T,  leading  to  considerable  translational 
freedom  for  the  ions.  This  should  result  into  quasi-lattices  with  at  best  only 
partial  long-range  order. 
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Table  2 - 1 . A c ompa  r i son  of  f^  obtained  from  Eqns . (2-14)  and  (2-15)  with 


experimental  values  for  NaCl 


Experimental 


0. 0005 


0.  932 


0.  948 


0.  975 


0.  001 


0.  915 


0.  935 


0.  965 


0.  002 


0.  894 


0.919 


0.  952 


This  first  approximation  to  the  partial  long-range  order  model  will  serve 
as  a basis  for  further  developments  shown  in  the  following  chapter  and  to  illus- 
trate concepts. 

In  this  model  we  will  retain  the  assumption,  later  to  be  removed,  of  ions 
occupying  on  the  average  the  centers  of  equal  cubes, but  will  consider  a partial 
degree  of  long-range  order  which  is  a function  of  temperature  and  concentration. 
The  temperature  dependence  is  analogous  to  that  of  binary  alloys  (e.  g.  Kittel, 
I959,  p.  340).  We  introduce  the  concentration  dependence  because  it  is  con- 
ceivable that,  as  the  ions  become  closer  together,  their  coulombic  energy  of 
interaction  increases  and  their  translational  freedom  decreases.  This  should 
lead  to  an  increase  in  the  degree  of  order  because  iF^/RTm  increases,  m is 
introduced  in  the  denominator  because  then  AE^m  is  the  energy  per  mole,  the 
quantity  to  be  compared  to  RT. 

The  concept  of  long-range  order  is  appealing  because,  on  the  average, 
the  nature  of  the  ions  which  surround  a given  central  ion  should  be  the  same  as 
that  of  the  ions  which  surround  the  ion  next  to  the  central  one,  and  so  on,  lead- 
ing to  a quasi-lattice.  In  Fig.  8-1  we  illustrate  this  process  for  the  perfect 
order  case. 

Let  us  now  define  a long-range  order  parameter  defined  by  each  pair  (T, 
m)  and  with  the  range  of  values 


0 < L ^ 1 


(3-1) 


3-2 


Fig.  3-1.  The  building-up  of  a quasi-lattice  in  two  dimensions  with  the  central 
ion  being  indicated  by  a circle. 
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so  that  Eqn.  (2-10)  becomes 


2„  1/3 

V-  AM°  L 

(1+0.  001  mM^) 


1/3 


4/3 

m 


(3-2) 


L = 0 for  T - ® or  m - 0 and  L = 1 for  T - 0 or  m - ®.  It  is  logical  that  AEj>  be 
zero  when  L = 0 (perfect  randomness)  because  then  there  is  an  equal  probability 
of  a cation  C or  an  anion  A occupying  any  given  lattice  site.  Any  one  ion  is, 
therefore,  surrounded  by  the  same  numbers  of  cations  and  anions  in  the  nearest 
neighbor  and  in  further  layers  so  that  the  net  potential  on  the  central  ion  is  zero. 

Let  us  next  express  L in  terms  of  the  degree  of  order  and  of  related  quan- 
tities. We  will  be  able,  by  means  of  an  artifice  and  of  the  Madelung  constant  to 
obtain  a rigorous  expression  for  AF^  in  contrast  to  the  case  on  non-ionic  binary 
alloys  for  which  successive  approximations  for  the  interactions  of  nearest 
neighbors,  next  nearest  neighbors,  and  so  on,  are  required. 

Let  us  consider  two  sub-lattices;  c for  the  cations  and  a for  the  anions, 
and  let  w be  the  fraction  of  'wrong  ions'  (C  on  a and  A on  c),  that  is,  the  fraction 
of  the  ions  of  a given  type  which  are  in  the  countcr-ion  lattice.  Let  us  set  w = 
0.25  for  example,  in  order  to  clarify  our  thoughts.  Then,  on  the  average,  one 
obtains  the  pattern  shown  in  Fig.  3-2.  This  figure  represents  the  average  state 
of  what  in  reality  is  a flickering  series  of  structures. 

Formally  we  can  then  represent  the  partially  ordered  lattice  by  a perfectly 
ordered  one  in  which  the  ionic  charge  is  (0.75  - 0.  25)ze  = (1  - 2 w)zc  and  one 


Fig.  3-2.  Average  lattice  configuration  for  the  fraction  of  ions  in  the  counter- 
ion lattice  being  w = 0.  25. 
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may  construct  Table  3-1.  d and£,  the  degrees  of  disorder  and  of  order,  are 
defined  by 


d = 2w  t = ( 1 - 2w) 


(3-3) 


Randomness  occurs  when  there  is  an  equal  probability  of  any  lattice  being 
occupied  by  an  anion  or  by  a cation  so  that  the  time-average  net  charge  at  any 
lattice  site  is  zero  (Fig.  3-3). 

Note  that,  if  the  effective  charge  is  (1  - 2w)ze  then  Eqn.  (2-10)  becomes 


2 1/3  2 

y amp  (i  - 2w)  4/J 

1/3  m 

(1  + 0.  001  mM2) 


(3-4) 


or 


(1+0.  001  mM2) 


1/3 


m 


4/3 


(3-5) 
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Tabic  3-1.  Values  of  w,  d = degree  of  disorder,  andf  - degree  of  order. 

w d / 


0.  5 

] 

0 

randomness 

0.4 

0.  8 

0.  2 

0.  3 

0.  6 

0.4 

0.  25 

0.  5 

0.  5 

0.  1 

0.  2 

0.  8 

0 

0 

1.0 

perfect  order 
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super-lattice  treatment.  There  are 


N m ( 1 - w)  cations  on  lattice  c or  anions  on  lattice  a 
A 


N m w cations  on  lattice  a or  anions  on  lattice  c 
A 


As,  from  Eqn.  (3-3) 


w = 0.  5 (1  -Z) 


1 - w = 0.  5 (1  +1) 


there  are 


0.  5 m ( 1 + £)  cations  on  lattice  c or  anions  on  lattice  a 


0 5 N m (1  -t)  cations  on  lattice  a or  anions  on  lattice  c 
A 


The  number  of  arrangements  of  these  ions  is 


N m! 
A 


. m y 

g(i>  = '[0.5  ( 1 + t)  N m]  I [0.  5 (1  -l)  N m]J 


From  Stirling' s theorem 


lnxl  = x(lnx  - 1 ) 


Eqn.  (3-8)  becomes 


(3-7) 


(3-8) 


(3-9) 


in  g(C)  ? 2N  m In  2 - NAm  [(1  +£)  In  (1  +£)  + (!  -Z)  In  (1  -£)]  (3-10) 
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In  g(  0)  - 2N  m In  2 


(3-11) 


where  g(o)  c orresponds  to  complete  randomness. 

Roltzman's  definition  of  entropy  applied  here  is 


I = Vn 


(3-12) 


as  g(£)  is  the  multiplicity  of  the  state  defined  by£  , i.  e.  , by  a pair  (T,  m)  at  a 
given  pressure  (the  atmospheric  one  in  our  case).  Therefore, 


= kfln  g (l)  - In  g(0)]  = -Rm  [l  + 1)  In  (1  + £)+(!  - l)  In  (1  -£)]  (3-13) 


and  the  configurational  free  energy  AF^  = AE^  - TAS^  is 


„ 2.  1/3  2 

hz  am°  1 

iFc = " m 

( 1 + 0.  001  mM  ) ' 


4/3 

m + 


+ RTm  [(1  + £)  In  ( 1 +£)  + (!  - l)  In  ( 1 -£)] 


(3-14) 


The  value  of  £ at  any  given  m which  yields  the  most  stable  configuration 

is  that  one  for  which  AF»  is  at  a minimum.  This  value,  indicated  by£  , can 

•e  mm 

be  obtained  graphically  or  by  setting  oAFj, /c£  = 0,  which  yields 


on  J/3 

_j i + t _ 2 ]lz  amp  m 

L n 1 -t  ~ RT  ( 1 + 0.  001  mM  ) 

min  ' 2 


(3-15) 


3-  13 


r 

A graphical  pr esentat  ion  of  AF^  versus  t is  shown  in  Fig.  3-4  and  a few 
selected  values  of  AF^,  AE^,  and  shown  in  Table  3-2.  It  can  be  seen 

that  there  is  a critical  concentration  (somewhere  between  m = 0.4  and  0.5 
according  to  our  numerical  results  not  shown  here)  below  which  AF^  is  a mini- 
mum only  at  £ =0,  that  is,  below  which  the  interactions  among  ions  are  weak 
enough  compared  to  the  thermal  energy  to  cause  the  long-range  order  to  break 
down.  This  is  analogous  to  the  critical  temperature  in  alloys,  or  to  a critical 
temperature  in  the  case  under  study,  above  which  there  is  no  long-range  order 
because  the  thermal  energy  is  much  larger  than  the  lattice  energy.  It  will  be 
shown  in  the  next  section,  however,  that  such  a critical  concentration  may  not 
exist  in  solution  when  the  translational  motions  of  the  ions  and  other  factors 
are  taken  into  consideration. 

It  is  interesting  to  observe  in  Fig.  3-4,  as  is  to  be  expected,  that£ 

min 

increases  with  increasing  m.  This  occurs  because,  as  the  distance  between 
nearest  neighbors  decreases,  the  coulombic  configurational  energy  increases 
in  magnitude  and  tends  to  induce  added  order. 

Let  us  next  compare  the  values  of  AF.  with  values  of  AF  based 

c(min)  (exp) 

upon  experiments.  AF  can  be  obtained  from 

(exp) 


. m 


AF  = 2RT  f in  f dm 

(exp)  Jo  ±(exp) 


(3-16) 


If  we  use  Eqn.  (1-27)  to  represent  log  f then 

±(exp) 


Table  3-2.  Examples  uf  Ihe  variation  of  AF^,  and  AS^,  versus  f. 

m t AF  AE  AS 

l f < 


0.  1 


0.  ) 
0.  2 
0.  3 
0.  4 
0.  5 
0.  6 
0.  7 
0.  8 
0.  9 


9. 5720x10° 
3. 8794 x 107 


-1.5260x10' 
-6.  1040 


0.4  0.1 

0.  2 
0.  3 
0.  4 
0.  5 
0.6 
0.  7 
0.  8 
0.  9 

0.  5 0.  1 

0.  2 
0.  3 
0.  4 
0.  5 
0.6 
0.  7 
0.  8 
0.  9 

1.0  0.1 

0.  2 
0.  3 
0.4 
0.  5 
0.  6 
0.  7 
0.  8 
0.  9 


8.  9248 
1. 6381xl08 
2.6708 
4.  06  29 
5. 93  1 2 
8.  4826 
1. 2!64xl09 
3. 1716x10^ 

1 . 4 7 1 0 x 1 0 7 
4. 0944 
9. 3365 
1. 9041xl08 
3.6096 
6. 5177 
1. 1456xl09 
2.  0211 
-4.  9906x1 06 

- 1. 743 2x1 07 
-2.  9415 

-2.  6574 
1. 4134 
1. 2882xl08 
3. 7592 
8.  5887 
1. 80 1 Ox  1 09 
-7.  309  I >d  08 
-2.  8730 
-6. 2682 
-1. 0629x109 

- 1 . 5495 
-2.  0143 
-2.  3405 
-2.  3213 

- 1. 5099 


-2.  4416 
-3. 8150 
-5. 4936 
-7.  4774 
-9. 7664 

- 1. 236  1 x 1 09 
-9*  6 156x  107 
-3.  8463x10s 
-8.6541 

- 1.  5385xl09 
-2.  4039 
-3.4616 
-4.7117 

-6. 1540 
-7. 7887 
-1. 2915x10s 
-5. 1660 
-1. I624xl09 
-2.  0664 
-3. 2288 
-4.6494 
-6. 3284 
-8.  2656 

- 1 . 046 1 x 1 0 1 0 
-3.  2141xl0R 

- 1. 2856X109 
-2.  8927 

-5.  1426 
-8. 0353 

- 1.  1 57 1 X 10  1 0 

- 1. 5749 
-2. 0570 
-2.6034 


-8.  328x  104 
-3. 3483X105 
-7. 5995 
-1. 3683X106 
-2.  1753 
-3. 2052 
-4. 4971 
-6.  1 205 
-8.  2252 
-3. 3314X105 
-1. 3393xl06 
-3. 0398 
-5.4731 
-8.  7011 
-1. 2821 x 1 07 

- 1. 7988 
-2.  4482 
-3. 2901 

-4.  1642x10S 
-1. 6742X106 
-3.7998 
-6. 84)4 

- 1 . 0876x  107 

- 1. 6026 
-2.  2486 
-3. 0603 

-4. 1 126X107 
-8. 3284xl0S 
-3. 3483x1 0^ 
-7. 5995 
- 1 . 3683X107 
-2.  1753 
-3. 2052 
-4.4971 
-6.  1 205 
-8.  2252 
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C 2 1/2 

Ll\.  = 2 x 2.  3026  RT  [ ^ {0.  5 B m - B m ' + 

(mterp)  3 C C 

C 


1/2  2 

+ 2.  3026  In  ( 1 + B „m  )}  + 0.5  C m + 


+ 0.4  Dm2’  3 + 0.333  E m3] 

v>  C 
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As  an  alternative,  AF,  , can  be  calculated  from 

(exp) 


AF(exp)  ” ^2°  + nNaCl  NaCl 


(3-18) 


where  n T _ and  n are  the  numbers  of  moles  in  the  given  amount  of  solution. 

1+2°  NaCl 

~ is  the  chemical  potential  of  the  water  and  App  is  the  configurational 

H^CJ  NaCl 

chemical  potential  of  NaCl.  As 


^NaCl=2RTlnWl 


(3-19) 


1 - 0 


- h2°  1000 
H O 2 m RT  _ 
11 2° 


(3-20) 


Eqn.  (3-18)  becomes 


3-18 


4IW55-5i(1  •’’h.o’ 


+ 2m  RT  In  f 


±NaCl 


(3-21) 


0 ^ and  M represent  the  practical  osmotic  coefficient  and  the  molecular 

H2  2 

weight  of  water.  Values  of  A F,.  , and  of  A F,  . are  compared  to  minimum 

b (interp)  (exp) 

values  obtained  from  Eqn.  (3-14)  in  Table  3-3.  It  is  obvious  that  the  theory  in 
its  present  form  does  not  work,  as  it  does  not  yield  correct  configurational 
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Tabic  3-3.  Comparison  of  AF  values  in  ergs  derived  from  experimental 

activity  coefficients  of  NaCl  at  25°  C with  minimum  values  of  AF^ 
obtained  from  Eqn.  (3-14).  See  Fig.  3-4. 


( inter p) 
Eqn.  (3-17) 


(exp) 

Eqn.  (3-21) 


AFf(min) 

Eqn.  (3-14) 


0.  001 

0.  005 

*7 

0.  01 

-3. 508x10 

0.  02 

-9. 53  2x 1 07 

0.  03 

-1. 700X108 

0.  04 

-2. 553xl08 

0.  05 

-3.493x10s 

0.  06 

-4. 505x 1 08 

0.  07 

-5.  580xlC8 

0.08 

-6. 771xl08 

0.  09 

-7. 891x10s 

0.  10 

-9. 1 16xl08 

0.  20 

-2. 318xl09 

0.  30 

-3. 945xl09 

0.40 

-5. 076xl09 

0.  50 

-7. 569X109 

0.  60 

-9.498X109 

0.  70 

-1.  148x 10 1 0 

0.  80 

-1. 350xl0]0 

0.  90 

-1. 556xl010 

1. 00 

-1.764xl010 

2.  00 

(*) 

Approximate  value  as  d 

-3. 516x10 


■1. 763x1 01U 
■3. 834x1 0 1 0 


-2.  942x10 
-1.  190xl08 
-5. 559x 1 08 
-1. 027xl09 
- 1. 65 1 x 1 09 
-2.  34  1 x 1 09 


0. 2-0. 4 
0.  2-0.4 
0.5-0. 7 
0. 6-0. 8 
0.6-0. 8 
0.6-0.  8 


Thus,  for  example,  found  a minimum  AF,  , at  1 = 0.  3 for  m = 0.  5 and 

(theor.  ) 

did  not  seek  a more  precise  value  by  adding  a second  decimal  figure  to  values 


of  t between  0.  2 and  0.  4. 
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4.  Partial  Lonn-^ai)Rc  Order  and  (lie  Parameter  v 
a)  The  parameter  v 

Several  important  factors  were  not  considered  when  Eqn.  (3-14)  was 
derived,  namely,  the  effect  of  ionic  motion,  possible  changes  in  dielectric 
constant  with  concentration,  repulsive  forces  caused  by  dielectric  constant 
gradients  (Dahe,  1972),  the  formation  of  ion-pairs  if  such  do  occur  in  chloride 
solutions,  and  the  change  in  effective  concentration  with  hydration. 

In  this  chapter  we  will  consider  the  effects  of  the  first  process  mentioned 
above  because  it  pertains  to  the  term  and  will  initially  lump  it  into  a 
parameter  v such  that  Eqn.  (3-14)  now  becomes 


w 2 1/3.2 

V -VP  1 * 4/3  , 

LF  = . — m + 

(1  + 0.  001  mM^) 


+ RTm  [(1  +1)  In  (1  + £)  + (1  -l)  In  (1  -£)]  (4-1) 

In  this  chapter  we  will  calculate  values  of  v and£  and  in  the  next  chapter 
we  will  recalculate  them  after  considering  hydration  corrections.  The  remaining 
processes  will  only  be  treated  indirectly  in  this  work.  Wc  will  proceed  in  this 
manner,  by  a step-wise  development  of  the  theory  for  two  reasons;  greater 
clarity  and  the  case  of  backtracking  if  in  the  future  one  of  our  steps  should  require 
revision.  After  recalculating  v and£,  we  will  attempt  to  justify  the  values  obtained 
by  considering  the  contributions  of  ionic  motion  and  of  dielectric  constant  effects. 


It  is  by  no  means  intuitively  clear  that  ionic  motions  around  central 
positions  (translations  in  dilute  solutions  tending  towards  vibrations  in  very 
concentrated  solutions),  and  we  will  justify  next  the  need  to  take  them  into 
consideration.  Imagine  if  you  will  the  three  ions  shown  in  Fig.  4-1,  with  the 
central  one  oscillating  between  its  two  nearest  neighbor  counter-ions. 

1/ r is  the  key  quantity  for  energies  of  interaction  and,  for  a given  x is, 
as  is  obvious  from  Fig.  3-5  t 


and 


with 


(4-2) 


r 

o 

r 

x 


1 

3 - (x/rQ)2 


(4-3) 


r > r 
o x 


(r  - b)2 
o 


r 

o 


(4-4) 


Thus, 


one  observes  that  AF^ 


is  increased  by  the  motion  of  ions  and  should  be 


"The  argument  that  follows  is  an  approximate  one  as  is  the  calculation  of  7 
in  section  4d  because  the  anions  are  assumed  fixed.  This  docs  not  detract 
fro?n  the  exactness  of  v calculated  from  experimental  data  but  only  from 
its  quantitative  intcrr>retatir»n 


Fig.  4-1.  Oscillating  ion.  The  broken  arc  sphere  represents  the  central 
cation  displaced  due  to  its  thermal  energy,  r is  the  average 
nearest  neighbor  distance,  b is  the  distance  of  closest  approach, 
and  r^  is  the  effective  (potential)  nearest  neighbor  distance  (see 
text),  x is  the  displacement  from  the  central  position. 


in  order  to  replace  r used  in  the  original  derivation  of  AE^,  (see  Eqn.  (2-2)). 
v in  Eqn,  (4-1)  is  now  replaced  by 


r , 

o 

V = =—  V 

r 

x 


(4-6) 


where  r is  the  actual  average  distance  belw'een  ions  and  v reflects  dielectric 
x 

effects.  7 will  be  calculated  in  a latter  chapter,  after  we  have  considered 
x 

the  effects  of  hydration.  For  the  time  being  we  will  limit  ourselves  to 

obtaining  values  of l and  v from  a comparison  of  Eqn.  (4-1)  with  AF(  v (see 

(exp) 

Eqn.  (3-16)  and  Table  3-3). 

b)  Procedure  for  the  calculation  of  t and  v 

The  procedure  for  the  calculation  of  -C  and  v was  as  follows: 

1.  At  each  value  of  the  molality  m the  values  of  AF  were  calculated  for  a 
series  of  values  of  t and  v by  means  of  Eqn.  (4-1). 

2.  For  each  (m,  v)  pair  there  was  a non-zero  minimum  value  of  AF^  provided 
that  v was  equal  to  or  larger  than  a critical  value  v^.  For  m = 0.  1 as  an 
example,  the  non-zero  minimum  value  of  AF^,  appeared  when  v ^ 4.0  and 
zero  was  the  only  minimum  value  for  v ^ 3.  0.  We  already  saw  in  the  pro- 


ceeding chapter  that  at  m = 0.  1 and  v 1.0  there  was  breakdown  of  long- 


4 -6 


range  order  and  now  we  see  that  Ihc  same  is  true  for  v - 2.  0 and  v - 3.0. 
Only  one  (m,  v)  pair  and,  consequently  one  value  of t yielded  a minimum, 
indicated  by  AF £(min)>  which  was  in  addition  equal  to  AF (exp)- 

O/ 

3.  We  found  the  values  of  AF^,  which  bracketed  AF^(min)  and  interpolated  the 

>|c 

corresponding  values  of  C and  v to  obtain  (,  and  v which  corresponded  to 

af  The  results  are  shown  in  Table  4-1  for  NaCl  and  in  Table  4-2  for 

£(min) 

KC1. 

As  is  to  be  expected,  the  absolute  value  of  AE^  increases  with  the 

concentration  because  the  ions  are  closer  together.  AS^  becomes  more 

negative  with  increasing  m because  it  reflects  the  change  of  entropy  in  going 

from  a random  configuration  to  a progressively  more  ordered  one.  The 

absolute  value  of  AF'\  . . = AF  , increases  with  concentration  and  so  does 

£(min)  (exp) 

the  degree  of  order  l" . The  parameter  v decreases  partly  due  to  the  decrease 

in  r /r  with  increasing  energy  of  interaction  which  tends  to  dampen  the  ionic 
o x 

motions.  This  can  be  seen  from  Eqn.  (4-3)  rewritten  in  the  form 


r 

o 

r 

x 


1 - (x/rQ) 


(4-7) 


As  the  energy  of  interaction  increases  due  to  the  increasing  concentration,  x/ 

r is  expected  to  decrease  from  1 at  infinite  dilution  (i.  e.  r / r - ®)  towards 
o ox 

q i r /7  ] ).  The  concentration  dependences  of  £ and  v are  plotted  in  figs. 

' o x 

4-2  and  4-3. 


C 


....... 


^ • — _ 
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Tabic  4-1.  Values  of  AF„  which  bracket  AF,  . corresponding  values  ofb 

i (exp) 

# '!•  >!< 

and  v,  and  interpolated  values  of  AF„,  . , = AF,  v,  -f  and  v . 

t(  mm)  (exp) 

The  quantities  corresponding  to  asterisks  are  shown  in  parentheses. 


The  results  arc  for  NaCl  at  25°  C. 


m 

AEf 

LSl 

LFt 

l 

V 

0.  005 

-9.067X107 
-9. 376x 1 07 
(-9. 366xl07) 

-2. 630X105 
-2.  712xl0S 
(-2. 709xl0S) 

-1. 224x1 07 
- 1. 288x 107 
(-1. 286xl07) 

0.  750 
0.  760 
(0. 760) 

5.72 
5.  76 
(5.76) 

0.  01 

-2.  172xl08 
-2.  176xl08 
(-2. 175xl08) 

-6. 121xl0S 
-6. 121xl0S 
(-6.  1 2 1 x 1 05 ) 

-3.47  1 X 1 0 7 
-3. 5l6xl07 
(-3. 508xl07) 

0.  800 
0.  800 
(0. 800) 

4.  78 
4.  79 
(4.79) 

0.  02 

-5. OOlxlO8 
-5.073X108 
(-5. 009X108) 

- 1. 358xl06 

- 1. 378xl0S 
(-1. 360X106) 

-9.5  17xl07 
-9. 642xl07 
(-9. 532xl07 

0.  835 
0.  840 
(0.  836) 

4.01 
4.  02 
(4.01) 

0.  03 

-8. 125xl08 
-8. 243x10s 
(-8.  212X108) 

-2.  1 6 1 x 1 06 
-2. 1 93x  1 06 
(-2.  1 85  x 1 0S 

-1.683xl08 
-1. 706xl08 
(-1. 700x10s 

0.  855 
0.  860 
(0. 859) 

3.62 

3.63 
(3. 63) 

0.  04 

-1. 15 2X1 09 
-1.  159X109 
(-1. 153X109) 

- 3 . 0 1 1 X 1 0 6 
-3. 01 1 X 1 06 
(-3. 01 1X106) 

-2.  547x10s 
-2.  6 1 5 x 1 0 8 
(-2. 553X108) 

0.  870 
0.  870 
(0.  870) 

3.  38 
3.40 
(3.  38) 

0.  05 

-1. 503xl09 
-1. 512xl09 
(-1. 505xl09) 

-3. 876X108 
-3. 876xl0S 
(-3. 876x10s) 

-3.470x108 
-3. 564x30s 
(-3.493xl08) 

0.  880 
0.  880 
(0. 880) 

3.  20 
3.  22 
(3.  20) 

0.  06 

-1. 865X109 
-1. 876xl09 
'-1. 874xl09) 

-4. 780X106 
-4. 790xl0S 
(-4. 789xlOS) 

-4.453x10s 
-4. 514xl0S 
(-4. 505x10s) 

0.  889 
0.  890 
(0. 890) 

3.06 

3. 07 
(3.07) 

-2.  2 1 Ox  1 09 
-2.  276xl09 
-2.  2 7 1 x 1 0 9 ) 

-5.  5S9X.106 
-5. 758xl0S 
(-5. 746xlOS) 

-5. 440x 1 08 
-5. 59 1x10s 
(-5. 580x10s) 

0.  890 
0.  900 
(0. 899) 

2.  94 
2.  96 
(2.96) 

ri.Sx  in9 
-~o,io9 
• x I09) 

-6.600xl0S 
-6. 600x10s 
(-6. 600xl0S) 

-6.694x10s 
-6. 740xl08 
(-6.717x10s) 

0.  90  1 
0.  901 
(0. 901) 

2.  86 
2.  87 
(2. 87) 
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m 

AE^, 

AS{ 

l 

V 

0.09 

-2.  987xl09 
-3. 075x1 09 
(-3. 026xl09) 

-7.406x10° 
-7.627x10° 
(-7. 503x10) 

-7. 795xl08 
-8. 013X108 
(-7. 89 1 X l O8) 

0.  900 
0.  910 
(0. 904) 

2.  78 
2.  80 
(2.79) 

0.  10 

-3. 437xl09 
-3.463x109 

(-3. 438x 109) 

-8.475x10° 

-8.475x10° 

(-8.475x10°) 

-9.  1 04x  1 08 
-9. 357X108 
(-9. 1 I6xl08) 

0.910 
0.  910 
(0. 910) 

2.  72 
2.  74 
(2.  72) 

0.  20 

-7. 630xl09 
-7. 696X109 
(-7.687xl09) 

- 1. 8 1 Ox  1 07 
-1. 810xl07 
(-1. 810X107) 

-2.  26  1 x 1 09 
-2.  327X109 
(-2. 318X109) 

0.  930 
0.  930 
(0. 930) 

2.  30 
2.  32 
(2.  32) 

0.  30 

-1. 219x10 10 
-1. 231X1010 

(-1. 229X1010) 

-2.  786xl07 
-2.  786 xlO7 
(-2.  786xl07) 

-3. 883X109 
-3. 999X109 
(-3. 945X109) 

0.  940 
0.  940 
(0.  940) 

2.  10 
2.  12 
(2.  12) 

0.40 

- 1 . 6 14x  1 0 1 0 
-1.623X1010 
(-1.622x10*°) 

-3. 714xl07 
-3. 7 14x 107 
(-3. 714xl07) 

-5. 069X109 
-5.  154xl09 
(-5. 076x1 09) 

0.  940 
0.  940 
(0. 940) 

1. 90 

1.91 
(1.91) 

0.  50 

-2.  177xl0]0 
-2.  210X1010 
(-2. 202xl010) 

-4.  822X107 
-4. 853X107 
(-4. 845xl07) 

-7. 395x1 09 
-7. 629X109 
(-7. 569xl09) 

0.  950 
0.  950 
(0. 950) 

1.  86 
1.  86 
(1. 86) 

0.  60 

-2. 639x10*° 
-2. 725xl010 
(-2.711x10*°) 

-5. 749x 1 07 
-5. 938x 107 
(-5. 907xl07) 

-9. 248x 1 09 
-9. 548xl09 
(-9.498X109) 

0.  950 
0.  960 
(0.  958) 

1. 78 
1. 80 
(1.80) 

0.  70 

-3.  1 90x  1 0 1 0 
-3. 227  x 1 0 1 ° 
(-3. 2 1 3x 1 0 1 °) 

-6. 927xl07 
-6. 927X107 
(-6. 927xl07) 

-1. 125x10*° 
-1. 162x10*° 
(-1. 148x10*°) 

0.  960 
0.  960 
(0. 960) 

1.72 
1. 74 
(1.73) 

0.  80 

-3.670X1010 
-3. 714x10*° 
(-3. 710xl010) 

-7. 917xl07 
-7. 9 1 7x 1 07 
( -7. 9 1 7x 1 07) 

-1.309x10*° 
-1.354x10*° 
(-1. 350X10*°) 

0.  960 
0.  960 
(0.960) 

1 . 66 
1.68 
(1.68) 

0.  90 

-4. 182xl010 
-4.  23 2x  1 0 10 
(-4.  21  2x  1 0 1 °) 

-8.  906x1 07 
-8.  906xl07 
(-8. 906xl07) 

-1.525x10*° 

- 1. 576x10*° 
(-1.556x10*°) 

0.  960 
0.  960 
(0. 960) 

1.62 

1.64 

(1.63) 

1.00 

-4. 680x10*° 
-4.  739X1010 
( - 4 . 7 1 1 X J 0 1 0 ) 

-9. 896x1 07 
-9. 896x107 
(-9. 89° X 1 07) 

- 1. 730x10*° 
-1.789x10*° 
(-1. 764x10*°) 

0.  960 
0.  960 
(0. 960) 

1.58 
1.60 
(1. 59) 

2.  00 

-o. 79 1 X 10  1 0 

- 1 . 004x10*  1 
862X10*°) 

-2.012X108 
-2.  046x1  08 
(-2.  022x  1 08) 

-3. 792x10*° 
-3. 940x10*° 
(-3. 834x10*°) 

0.  970 

0.  970 
(0.  970) 

1. 33 

1. 35 
(1.34) 

Tabic  4-2.  Counterpart  of  Table  4-1  but  for  KC1  solutions  at  25'  C. 


m 

AFl 

i 

0.  0001 

-3. 609x1 04 
-4 . 0 1 9x 1 04 
(-3. 850xl04) 

0.  510 
0.  520 
(0. 516) 

0.  0005 

-4.012xl05 
-5. 028xl05 
(-4. 253xl05) 

0.  600 
0.  640 
(0.  609) 

0.  001 

-1. 020xl06 
-1. 302X106 
(-1.  193xl06) 

0.640 
0.  670 
(0.  658) 

0.  005 

-1. 191X107 
- 1 . 353  X 1 07 

(-1. 286  X 1 07 ) 

0.  750 
0.  770 
(0.  762) 

0.  01 

-3. 517x1 07 
-3. 562X107 
(-3. 543x1 07) 

0.  800 
0.  805 
(0.  803) 

0.  02 

-9. 653x 1 07 
-9. 779  X 1 07 
( -9. 667x 107) 

0.  839 
0.  840 
(0. 839) 

0.  03 

-1. 709xl08 
- 1 . 73 1 Xl08 
(-1. 729X108) 

0.  859 
0.  860 
(0. 860) 

0.  04 

-2.  553X108 
-2.622x10s 
(-2. 605x10s) 

0.  871 
0.  873 
(0.  873) 

0.  05 

-3.  5 28x 1 08 
-3. 575x10s 
(-3. 573x10s) 

0.  885 
0.  885 
(0. 885) 

0.  06 

-4. 592X108 
-4.654x10s 
(-4.620xl0S) 

0.  891 
0.  893 
(0.  892) 

17.  900 

18. 000 
(17.  959) 

11. 000 
1 1. 200 
(1  1. 047) 

8.  900 
9.  100 
(9.  023) 

5.  700 
5.  800 
(5.  759) 

4.  790 
4.  800 
(4.796) 

4.  020 
4.  030 
(4.  021) 

3.  630 
3.  640 
(3.639) 

3.  380 
3.400 
(3.395) 

3.  210 
3.  220 
(3.  220) 


3. 0  80 
3. 090 
(3. 085) 
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m 

' l 

i 

V 

0.  07 

-5.692X108 

0.  899 

2.  980 

-5. 769x10* 
(-5.734x10  ) 

0.  900 

2.  990 

(0.  899) 

(2. 985) 

0.  08 

-6. 865xl08 

0.  903 

2.  880 

-6.  958xl08 

0.  905 

2.  890 

(-6. 9 1 Ox  1 08) 

(0. 904) 

(2. 885) 

0.  09 

-8. 056xl08 

0.  907 

2.  800 

-8.  1 66  X 1 0 8 

0.  909 

2.  810 

(-8.  1 4 1 x 1 08) 

(0. 908) 

(2. 805) 

0.  10 

-9.4] 1 x 108 

0.  915 

2.  740 

-9. 5 3 9 X 1 0 8 

0.  915 

2.  750 

(-9.422X108) 

(0. 915) 

(2.  741) 

0.  20 

-2.419X109 

0.  935 

2.  340 

-2. 453X109 

0.  936 

2.  350 

(-2.434x109) 

(0. 935) 

(2.  344) 

0.  30 

-4. 178xl09 

0.  947 

2.  140 

-4. 237 xlO9 

0.  9 48 

2.  150 

0.  40 

(-4.  1 94  x 1 09) 
(-6.  138xl09) 

(0. 947) 

(2. 143) 

0.  50 

-8. I65xl09 

0.  960 

1.910 

-8. 285xl09 

0.  960 

1. 920 

(-8. 221 XlO9) 

(0. 960) 

(1.915) 

0.  60 

-1. 043xl010 

0.  962 

1. 840 

(- J . 042x10 10) 

(0.  962) 

( 3. 840) 

0.70 

-1. 255xl010 

0.  964 

1. 770 

-1. 274xl010 

0.  965 

1 . 780 

( - 3 . 27  1 X 1 0 1 °) 

(0. 965) 

(1.778) 

0.  80 

- 1 . 492x 1 0 1 0 

0.  966 

1.  720 

- 1.  537xl010 

0.  968 

1. 740 

(-1. 507xl010) 

(0. 966) 

(1.721) 

4-  n 


m 

AF> 

l 

V 

0.  90 

- 1 . 7 1 8x 1 0 1 0 

0.  967 

1.670 

- 1 . 77 1 x 1 0 1 0 

0.  971 

1. 690 

( - 1 . 75 ] x 1 0 1 °) 

(0. 970) 

(1.686) 

1.0 

- 1 . 866x 10 1 0 

0.  965 

1. 600 

-2.019X1010 

0.  970 

1.650 

( - 1 . 992x 1 0 1 °) 

(0. 969) 

(-1.641) 

2.  0 

-4. 642x 1 0 1 0 

0.  978 

1.400 

-4. 798x1010 

0.  980 

1. 420 

(-4.657X1010) 

(0.  978) 

(1.402) 
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One  may  next  fit  the  curves  in  Figs.  (4-2)  and  (4-3)  with  empirical 
equations  which  are  then  introduced  into  (4-1).  An  equation  for  log  f can 
then  be  obtained  from  the  relation 


$ 

log  f = ! 5JLQ ™ \nl 

h ± 2 x 2.  3026RT  Bm 


(4-8) 


We  did  not  attempt  a fit  over  the  whole  range  of  concentrations  as  the 
main  purpose  of  this  work  so  far  is  to  illustrate  concepts  and  to  provide  a 
model  rather  than  to  predict  activity  coefficients.  We  found  that  for  m > 0.  03 


0.9683m 
(0.  00580  + m) 


(4-9) 


provided  an  excellent  fit  with  an  R = 0.  99999.  For  m > 0.  1 


v = 1.592/m 


0. 2356 


(4-10) 


had  an  R = 0.  9999.  The  equation  for  log  f then  became 


l0g  f±  2.  3026  RT 


, * 0. 00 1M 

, rrJL  2.  3 56  0.  183  2 21  2 

3m  " m " 3p  m " 0.  9682m  " 3(  1 + 0.  00  ImM. 


]]- 


T&S 


m 


+ RTm  [In  (1  + C)  - In  ( 1 -C)]  - ~~~ 


*2 


m 0.  9682m 


(4-11) 


Resulting  values  of  f are  compared  to  experimental  ones  in  Table  4-3. 
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Tabic 

4-3.  A comparison  of  f obta 

ined  from  Eqn. 

(4-11)  with  cx 

peri  mental 

values 

f for  NaCl 

i(exp) 

at  25°  C 

m 

f 

± 

f±(cxp) 

m 

f±(  exp) 

0.  1 

0.  776 

0.  778 

0.  6 

0.  674 

0.673 

0.  2 

0.  734 

0.  735 

0.  7 

0.667 

0.667 

0.  3 

0.  710 

0.  710 

0.  8 

0.662 

0.  662 

0.  4 

0.693 

0.  693 

0.9 

0.657 

0.  659 

0.  5 

0.683 

0.681 

1.0 

0.653 

0.657 
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A similar  procedure  can  of  course  be  usrd  for  m < 0.  1 clown  to  any  values  by 
finding  the  propor  fitting  equal  ions  fori  and  v . Table  4-3  is  presented  for 

*J;  5jc 

the  sake  of  completeness  but  does  not  prove  anything  as  we  oblained-C  and  v 

by  a curve  fitting  procedure  in  the  first  place. 

Rather,  our  contention  is  that  the  overall  logical  framework  of  the  model 

1 1 * * 

is  superior  to  that  of  the  Debye-Huckel  theory  and  that£  and  v follow  expected 
trends  and  fall  within  expected  ranges,  thus  providing  credibility  for  the  model. 
In  the  approach  to  infinite  dilution  we  would  of  course  obtain  the  same 

tl 

numerical  results  for  f as  those  provided  by  the  Debye-Huckel.  This  would 
occur  becaxjse  our  starting  point  to  obtain  f at  high  dilutions  would  be  the 
empirical  equation  (1-27)  which  assumes  the  limiting  form 


log  f 


±{  inter  p) 


A I 


0.  5 


(4-12) 


with  A = -0.  5108  (see  Table  1-3).  However,  in  our  case  the  equation  corres- 
ponding to  Eqn.  (4-11)  would  be  a more  complicated  function  of  m than  Eqn. 

. a. 

*i' 

(4-12)  and  would  depend  upon  the  equations  used  to  fit  empirically  l and  v as 

functions  of  m.  Note  that  Eqn.  (4-12)  was  obtained  under  the  false  premise 

that  n.  decays  to  the  bulk  value  n.  at  large  distances  from  the  central  ion  J 

(see  Eqn.  (1-2)  and  its  critique  in  section  6-e). 

The  reason  for  the  use  of  an  empirical  equation  for  l is  that  an  explicit 

r ieorous  value  of  versus  m cannot  be  obtained.  Rather,  in  analogy  to 

b min 

Eqn.  (3-15),  one  obtains,  by  setting  = 0, 
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, >♦*  2 V2AMp1/3„,'/3 

1 ^ rmn i M 

Vmin  ^min  ^ ^min  RT  (1  + 0.00]  lnM^)  ^ 


(4-13) 


The  subscript  "min"  refers  to  values  of  l and  v which  yield  minimum  values  of 

AF  which  do  not  necessarily  equal  AF  . If  they  do  then  the  symbols 

£(min)  (exp) 

# it- 

l , v and  A ^ are  used.  Furthermore,  the  symbols  expj  anc*  ^(interp) 

are  equivalent  as  Eqn.  (1-27)  is  used  to  obtain  the  values  of  f . = f . . 

xv  i nv  erp)  xv  cx  p j 

which  are  integrated  to  get  AF,  . It  is  straightforward  to  show  that  the 

( exp) 

i',Z  % 

values  of  Z and  v presented  in  parentheses  in  Tables  4-1  and  4-2  obey  Eqn . 
(4-13)  quite  closely  with  slight  deviations  due  to  the  interpolation  procedure 
used  in  the  Tables.  Thus,  at  m = 0.01  the  left  term  of  Eqn.  (4-13)  is  0.  573 
while  the  right  one  is  0.  572  and  at  m = 1 the  values  are  2.  55  and  2.  58  respec- 
tively. 

Another  reason  for  seeking  empirical  fits  at  this  point  is  that  the  term 
/ 

r /F  in  v = (r  /r  ) v is  quite  complicated,  as  will  be  shown  in  the  next 
ox  ox 

/ 

section,  and  v is  not  yet  known  as  a function  of  m. 

d)  Determination  of  r /F  and  v^  ^ 

o x 

The  concept  of  potential  distance  r^  was  introduced  in  section  4-a.  It 
is  the  interionic  nearest  neighbor  distance  to  be  used  in  the  eoulombic  energy  cf 
interaction  when  an  ion  is  displaced  by  x from  its  equilibrium  position,  at 
which  the  nearest  neighbor  distance  is  r . 

The  average  nearest  neighbor  distance  r^  resulting  from  the  balance 
between  the  energy  of  interaction  which,  combined  with  the  energy  of  hydration, 


L 


tends  to  maintain  a regular  distance  between  ions  and  the  thermal  energy  which 


• See  footnote  on  page  4-2. 
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produces  translation  (or  vibration  at  very  high  concentrations)  is  given  by 


r 

x 


r° 

| r exp  (-AEn/RTm)  dr 
I x <-  > 


min 


exp  (-iE,/RTm)  dr 
Jt  x 


(4-14) 


where  (see  Eqn.  (4-4)) 


r = r ( 1 
min  o 


(4-15) 


The  exponential  term  can  be  written  in  the  form 


AEp  . r 
-c(min)  o 


1 


RT  m 


^-=Kv  ^ 

r r 

x x 


(4-16) 


where 


K 


AE  ’ . r 

-c(min)  o 


(4-17) 


v 


x 


RT  m 


f 


AD-A031  316  OREGON  STATE  UNIV  CORVALLIS  SCHOOL  OF  OCEANOGRAPHY 

A DYNAMIC  PARTIAL  ORDER  MODEL  OF  ELECTROLYTE  SOLUTIONS. ( U) 

OCT  76  R M PYTKOWICZ*  K JOHNSON*  C CURTIS  N00014-76-C-0067 

UNCLASSIFIED  REF-76-13  NL 

• 


2 of  2 
A0A03I3I6 


END 


DATE 

FILMED 

12  -76 


4-21 


is  a constant  at  each  m.  The  values  of  AE„,  . , and  v can  be  obtained  for 

min) 

NaCl  at  25 c C from  Table  4-1.  The  distance  of  closest  approach  1)  in  Eqn. 
(4-15)  can  be  selected  as  a reasonable  value  for  hydrated  ions.  The  term  m 
appears  in  the  denominator  of  Eqn.  (4-17)  to  convert  AE^  to  a molal  basis  as 
RT  refers  to  the  thermal  energy  per  mole. 

Note  that  in  Eqn.  (4-14)  we  are  applying  the  Boltzman  distribution  to  a 


series  of  systems,  each  with  a given  energy  of  interaction  AE^ 


which  is 


a function  of  m,  rather  than  to  individual  interacting  particles  as  was  done  by 

It 

Debye-Huckel  and  successors. 

By  setting  y = 1/r  , Eqn.  (4-14)  becomes 


which  integrates  to 


| y exp  (K  v y)  dy 

1/r 


/ y exp  (K  v y)  dy 

1/r  . 

min 


(4-18) 


Kvy  ' Ky  ...  ',2 

{.s — - _ K^o_  + (Kxj_  fln  t 5 

,2  2y  2 .. 


' i 1/r 

? (Kvy)  ,, 

. . ixi!  )S  1/r 
1=1  min 


gK v y • « <k  v v)  i 

+Kv  fin  y + I .-"/HI  1/r 


' .1  ,/ro 


(4-19) 


fc. 
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/ 

The  problem  consists  in  finding  at  each  molality  m a value  of  v such  that 

/ 

r calculated  from  (4-19)  and  v obey  the  relationship 


(4-20) 


where  v for  NaCl  is  found  in  Table  4-1.  This  value  for  v is  not  unity  because 
/ 

for  v =1,  Eqn.  (4-20)  yields  the  data  shown  in  Table  4-4. 

/ 

We  will  not  seek  the  correct  values  of  v and  r at  this  point  because 

x 

/ 

first  we  wish  to  apply  hydration  corrections  in  the  next  chapter  to  render  v , 
r and  the  general  model  more  realistic. 


f 

Table  4-4.  Values  of  r , r from  Eqn.  (4-19)  with  v = 1,  (r  /F  ),  and 

ox  ox 

v from  Table  4-1  at  different  molalities  and  distances  of  closest 

o 

approach  b for  NaCl  at  25°  C.  Values  of  r and  b .n  A. 

_ * 1 


m 

r 

o 

r 

X 

(r  /r  ) 

O X 

V 

b 

0.  1 

20.31 

11. 07 

1.  83 

2.  72 

2.  76 

0.  1 

20.  3 1 

12.69 

1.  50 

2.  72 

3.  80 

1.0 

9.  479 

7.  73 

I.  23 

1.60 

3.  80 

5.  The  Hydration  Correction 


In  this  chapter  we  will  derive  a hydration  correction  in  a manner  similar 
to  that  derived  in  Robinson  and  Stokes  (1958),  with  a small  change  as  we  will 
use  molal  quantities  throughout,  and  will  apply  the  correction  to  our  lattice 
theory. 

a)  Derivation  of  the  hydration  equation 

The  free  energy  F of  a given  amount  of  solution  with  n^  moles  of  water 
and  one  mole  of  solute  which  dissociates  into  v moles  of  cations  and  moles 
of  anions  can  be  expressed  in  two  equivalent  ways: 


F = n p + V p + v p 
w w C C A A 


F = (n  - h)p  + V p + V p 
w w C C A A 


(5-2) 


where  h is  the  hydration  number  (number  of  moles  of  water  which  solvate  one 

/ 

mole  of  salt)  and  p.  is  the  chemical  potential  of  the  hydrated  form  of  ion  i. 
Eqns.  (5-  1)  and  (5-2)  are  identical  because  F is  independent  of  where  the 
water  of  hydration  is  assigned.  Thus;  setting  these  equations  equal  to  each 
other  and  developing  the  chemical  potentials  by  means  of 


p = p.  + RT  In  f.  m. 
ri  *1  ii 


(5-3) 


. ..  »\  I 





one  obtains  the  expression 


As  m = 1000/  18  n and  m = 1000/18  (n  - h) 


if  one  remembers  that  there  are  n moles  of  water  per  mole  of  salt  in  the 


solution 


Introducing  Ecjn.  (5-5)  into  (5-4)  one  obtains 


I u'wwei-v  " ■ 


5-3 


with  v = V + V . As  n — <*>  all  activity  coefficients  and  a tend  to  unity. 
C A w ’ w 

their  logarithms  lend  to  zero,  and 


n - h 
, w 

V In V In 


1 - h/n 


w 


1 


- 0 


(5-7) 


Therefore,  all  the  terms  with  standard  chemical  potentials,  which  do  not 


depend  upon  n , must  add  up  to  zero  and  Eqn.  (5-6)  is  reduced  to 


V In  f - v„  In  f + v In  f - v.  in  fA 
CC  CCAA  A A 


n - h 

r y In  + h In  a 

n w 


(5-8) 


w 


Then,  remembering  that 


VC  VA  1/V 
±=  (fC  fA  > 


(5-9) 


one  obtains 


Inf  = In  f + — In  a + In  (1  - 0.  018  h m) 
± ± V w ' 


(5-10) 





— 


t 
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as  n = 1000/18  m leads  to  (n  - h)/n  = (1  - 0.018  h m). 
w w w 

In  terms  of  the  practical  osmotic  coefficient 


, 5000  • 

0 = -77 In  a 

18  v m w 


(5-11) 


Eqn.  (6-10)  becomes 


Inf  = In  f - 0.018  h m 0 + In  (1  - 0.018  h m) 
i-  i 


When  one  deals  with  mole  fractions 


N_  = 


2 n + V 
w 


N 


2 n - h + V 

w 


and  the  counterpart  of  Eqn.  (5-5)  is 


N n - h + V 
2_  _ _w 

M ' ' 11  + V 

N w 

2 


(5-12) 


(5-13) 


(5-14) 


Hence  the  difference  between  our  results  and  those  presented  in  Robinson  and 
Stokes  (1958). 

jJ: 

The  activity  coefficient  that  corresponds  to  in  our  theory  is  f^ 

in  Eqn.  (5-12)  because  we  did  not  make  an  allowance  for  the  change  in  effective 
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m resulting  from  the  removal  of  water  of  hydration  in  going  from  infinite 

dilution  to  molality  m.  The  experimental  mean  activity  coefficient  does 

/ 

include  this  effect  and  corresponds,  therefore,  to  f . 

b)  Procedure  for  the  calculation  of  l and  v 

The  counterpart  of  the  procedure  in  section  4-b  is  now  as  follows: 

Firstly  one  calculates  f , the  mean  activity  coefficient  corrected  for 

±(corr) 

hydration,  by  means  of  Eqn.  (5-12)  in  the  form 


Inf,  = In  f + 0.018  h m 0 - In  (1  - 0.018  h m)  (5-15) 

i(corr)  ±(interp) 


where  f . , is  obtained  from  Eqn.  (1-27). 

£(  inter p) 

For  m < 0.  1,  f can  be  obtained  from  the  empirical  equation 

Ifinterp) 


A.  I'*5 

log  f = ° ^ - 

& + \ 0 5 

± 1 


(5-16) 


where  = -0.  5108.  In  the  case  of  NaCl  we  saw  that  the  best  fit  between  Eqn. 
(5-16)  and  the  experimental  data  is  obtained  for  = 4.60A.  As  = 0.3287 
(see  Table  1-1),  Bv  = BQ  * aQ  = 1.512. 

In  the  range  in  which  an  equation  of  the  type  (5-16)  is  valid,  the  practical 
osmotic  coefficient  can  be  calculated  from  (Robinson  and  Stokes,  )958,  p.  34) 


0=1- 


2.  3026  A^  m 


(5-17) 
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with 


o(x)  = ~ [d  + x)  - 2 In  (]  + x)  - ] 

x 


(5-18) 


and 


N 0.5 
x = B m 


(5-19) 


At  values  of  m i 0.  1 there  are  tabulations  of  0 versus  m. 

One  then  fits  f by  means  of  an  empirical  equation  and  calculates 

i(corr) 


AF 


(corr) 


m 

= 2RT  I In  f 


±(corr) 


d m 


(5-20) 


Values  of  AF„  as  functions  of  v and  t at  each  molality  are  obtained  as 

•t(min) 

before  either  by  using  printouts  of  AF^  given  by  Eqn.  (4-1)  for  arbitrarily 

selected  (£,  v)  pairs  or  by  selecting  values  of  v arbitrarily  and  calculating  l 

from  Eqn.  (4-13).  One  then  selects  that  AF„,  . labeled  AF  . which 
' l(  mm)  ^(m  in) 

equals  AF,  , and  thus  obtains  the  values  of  l and  v . 

(corr) 

In  contrast  to  chapter  4 in  which  there  was  only  one  pair  (c  , v ) for  each 
molality  m,  there  will  now  be  one  (t  , v ) pair  for  each  (m,  h)  pair,  where  h 
is  the  hydration  number,  that  is,  the  number  of  moles  of  water  which  solvate 


each  mole  of  salt.  Our  method  will  only  tell  us  which  value  of  h is  the  correct 


one  when  we  manage  to  derive  the  term  v (see  Eqn.  (4-6))  from  first  principles 
and,  at  present  wc  must  select  values  of  h based  on  literature  data, 
c)  Calculation  off  and  v 

We  will  illustrate  this  procedure  for  NaCl  solutions  at  25°  C.  The  first 

step  is  to  obtain  values  of  the  osmotic  coefficient  0 by  means  of  Eqn.  (5-17) 

for  m < 0.  09  and  from  the  literature  'Robinson  and  Stokes,  1958)  for  higher 

\ 

concentrations.  The  values  of  A and  B are  -0.5108  and  1.512  as  was  shown 

in  the  previous  section.  The  resulting  values  of  0 are  presented  in  Table  5-1. 

Next,  by  means  of  Eqn.  (5-15)  we  calculated  the  values  of  f as  a function 

i(corr ) 

of  m and  of  assumed  values  of  h between  1 and  10.  The  results  are  shown  in 
Table  5-2. 

Tables  5-3  and  5-4  contain  the  values  of  0 used  to  calculate  f for 

±(corr) 

KC1  and  the  resulting  values  of  f , Tables  5-5  and  5-6  contain  similar 

±(  corr) 

results  for  LiCl. 

The  activity  coefficients  in  Table  5-2  can  be  expressed  by  equations  such 
as  (1-27).  The  constants  from  such  equations  determined  at  each  value  of  h 
were  then  used  to  calculate  the  free  energy  corrected  for  hydration  A F 

(corr) 

Values  of  t , v and  AF,,  were  then  calculated  as  discussed  previously  for 

f.(min)  r 7 

each  m and  h.  A summary  of  the  results  for  NaCl  is  shown  in  Table  5-7. 

Tables  5-8  and  5-9  contain  the  results  for  KC1  and  LiCl. 


M 
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Tabic  5-1. 

Values  of 

the  practical 

osmotic  coefficient  0 for 

NaCl  at  25°  C 

versus  th< 

e molality  rn. 

m 

0 

m 

0 

m 

0 

0.  00(5 1 

0. 9962 

0.  05 

0. 9437 

0.  40 

0. 9203 

0. 0005 

0. 9917 

0.  06 

0. 9406 

0.  50 

0.  9209 

0.  001 

0.  9884 

0.  07 

0. 9379 

0.  60 

0.  9230 

0.  005 

0.  9762 

0.  08 

0. 9356 

0.  70 

0. 9257 

0.  01 

0.  9683 

0.  09 

0. 9336 

0.  80 

0. 9288 

0.  02 

0. 9587 

0.  10 

0.  9324 

0.  90 

0.  9320 

0.  03 

0.  9523 

0.  20 

0. 9245 

1.0 

0. 9355 

0.  04 

0.  9475 

0.  30 

0. 9215 

2.  0 

0. 9833 

Table  5-2.  Values  of  f , the  mean  activity  coefficient  corrected  for  hydration,  as  a function  of  m 
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Table  5-3.  Values  of  the  practic  al  osmotic  coefficient  0 for  KC1  at  25°  C 


versus  the  molality  m. 

m 

0 

m 

0 

m 

0 

0. 0001 

0.  9962 

0.  05 

0.  9407 

0.40 

0. 9017 

0. 0005 

0. 9916 

0.  06 

0.  9372 

0.  50 

0. 8989 

0.  001 

0. 9883 

0.  07 

0.  9342 

0.  60 

0. 8976 

0.  005 

0.  9758 

0.  08 

0.  9315 

0.  70 

0. 8970 

0.  01 

0.  9675 

0.  09 

0.  9292 

0.  80 

0. 8970 

0.  02 

0. 9572 

0.  10 

0.  9266 

0.  90 

0. 8971 

O 

o 

LO 

0.  9502 

0.  20 

0.  9130 

1. 00 

0. 8974 

0.  04 

0. 9450 

0.  30 

0.  9063 

2.  00 

0. 9124 

UP 

O' 

X 

o 

o 

IP 

CP 

X 

IP 

rj 

o 

ip 

X 

o 

PJ 

X 

CO 

e 

Pi 

*T 

x> 

X 

in 

ro 

ro 

X 

X 

QC 

X 

CP 

cr 

o 

PJ 

-* 

*— I 

PJ 

o 

CP 

X 

ro 

Pi 

O 

ro 

Pi 

f'- 

CO 

Pi 

X 

p 

X 

ro 

o 

X 

CP 

r— « 

ro 

p 

PJ 

X 

UP 

ro 

p* 

X 

X 

X 

CP 

CP 

O 

»— i 

-* 

~ 

o 

CO 

p] 

o 

X 

p 

o 

ro 

pj 

p 

PJ 

p 

p 

X 

X 

X 

o 

rO 

p 

pj 

p 

ro 

p- 

p 

X 

X 

• 

X 

CP 

CP 

O 

X) 

X 

p 

rH 

o 

X 

rf 

CP 

ro 

. — < 

r—t 

CO 

o 

*— « 

X 

CP 

X 

rr 

up 

X 

p- 

o 

rO 

X 

o 

UP 

p- 

p 

P 

• 

p 

X 

X 

X 

CP 

CP 

pi 

up 

CO 

o 

P- 

o 

O 

ip 

o 

p 

UP 

UP 

CP 

X 

X 

CO 

Pi 

CP 

fO 

*r 

X 

X 

r— < 

rf 

X 

P- 

p 

p 

p 

p- 

X 

X 

X 

CP 

X 

CP 

rO 

CP 

X 

o 

CP 

CO 

X 

Pi 

CP 

X 

•— « 

p 

o 

o 

»"  < 

cO 

Tf 

X 

CP 

< — ' 

p* 

p- 

p 

P 

p 

p- 

p 

X 

PJ 

ro 

p 

Pi 

p 

p- 

o 

ro 

O 

o 

ro 

r-i 

•—< 

o 

p 

CP 

CP 

CP 

<P 

o 

Pi 

*r 

UP 

°t 

x 

X/ 

X 

p 

• 

p 

• 

p 

p 

P- 

iP 

o 

P- 

X 

X 

X 

CP 

X 

CO 

cO 

— 1 

CP 

CP 

•— « 

X 

CO 

pi 

ro 

p- 

p 

X 

X 

p 

p 

X 

CP 

O 

X 

X 

X 

X 

X 

x 

X 

X 

P 

co 

LD 

o 

x 

o 

CP 

P) 

CP 

X 

** 

UP 

CP 

** 

IT- 

IP 

■*r 

UP 

x 

X 

X 

X) 

X 

X 

X 

X 

x 

P 

IP 

P 

CO 

X 

CP 

rr 

P- 

ro 

o 

CP 

a- 

rr 

cO 

Pi 

pj 

•— i 

*— « 

•— 1 

o 

o 

X 

o 

X 

X 

X 

X 

X 

X 

*t 

X 

O 

IP 

X 

't 

X 

X 

P 

O 

(O 

PJ 

p 

pi 

X 

rJ 

r—4 

o 

CP 

X 

p- 

p 

X 

X) 

X 

X 

• 

UP 

IP 

UP 

IP 

• 

up 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

I*- 

X 

a 

o 

PJ 

X 

X 

o 

r-H 

*— 1 

PJ 

Tabic  5-5.  Values  of  the  practical  osmotic  coefficient  p or  LiCl  at  25°  C 
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Tabic  5-7. 

The  values 

of  £ , v and  AF 
l 

for  NaCl  at  a 

range  of  molalities 

and  various 

h. 

m 

h 

Jjt 

l 

❖ 

V 

.0001 

0 

. 5203 

17.  95 

-3. 840x10^ 

2 

. 5203 

17.  95 

-3. 840x10* 

4 

. 5203 

17.  95 

-3. 840x10* 

6 

. 5203 

17.  95 

-3. 840x10* 

8 

. 5203 

17.  95 

-3. 840x10* 

10 

. 5203 

17.95 

-3. 840X104 

. 001 

0 

.6585 

9.  020 

-1. 190X10^ 

2 

. 6580 

9.  017 

- 1.  185x10? 

4 

.6580 

9.  017 

- 1.  185x10? 

6 

. 6580 

9.  017 

- 1.  185x10? 

8 

.6575 

9.  014 

-1.  180x10? 

10 

.6575 

9.  014 

-1.  180X1 06 

.01 

0 

. 7990 

4.  786 

-3.  504x10^ 

2 

. 7985 

4.  783 

-3.490x10^ 

4 

. 7978 

4.  779 

-3.  469x1 0^ 

6 

. 7970 

4.  774 

-3.449x10.1 

8 

.7965 

4.  771 

-3. 435x10.: 

10 

.7960 

4.  768 

-3.4  22x 1 0? 

. 1 

0 

. 9093 

2.715 

-9.  110X10® 

2 

. 9073 

2.  701 

-8.  929*10* 
-8. 752x10* 
-8. 580x10* 
-8.413x10* 
-8.  251x10 

4 

. 9053 

2.  837 

6 

. 9033 

2.673 

8 

.9013 

2.660 

10 

. 8993 

2.653 

. 5 

0 

. 9530 

1.  858 

-7.  582x10^ 

2 

.9485 

1. 820 

-7.  1 33x 1 0g 

4 

. 9435 

1. 782 

-6. 686x1 09 

6 

. 9383 

1.745 

-6.  266x10^ 

8 

. 9320 

1. 706 

-5.  820x1 0g 

10 

. 9250 

1.667 

-5.  378xl0V 

.0 

0 

. 9633 

1. 562 

-1. 766x10 
-1.  594X101* 
-1.417x10 

- 1.  238x10^ 

- 1. 058x10, ° 

2 

. 9565 

1.  505 

4 

. 94  80 

1.446 

6 

. 9373 

1.  386 

8 

.9235 

1.321 

10 

. 9050 

1. 253 

-8.  730x10 

it  Milifltaiilitoi  ■ 
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Table  5-H.  The  values  of  i , v and  for  KC1  at  a range  of  molalities 


and  various  h. 


m 

h 

l 

V 

0. 0001 

0 

. 5205 

17.  95 

-3. 848x10 

2 

. 5 205 

17.  95 

-3. 848x10 

4 

. 5205 

17.  95 

-3. 848x10 

6 

.5205 

17.  95 

-3. 848x10 

8 

. 5205 

17.  95 

-3. 848x10 

10 

. 5205 

17.  95 

-3. 848x10 

0.  001 

0 

.6585 

9.  022 

-1.  192x10 

2 

.6585 

9.  021 

-1.  190x10 

4 

.6585 

9.  021 

-1.  190x10 

6 

.6580 

9.  018 

-1.  185x10 

8 

. 6580 

9.  018 

-1.  185x10 

10 

.6580 

9.0)8 

-1.  185x10 

0.  01 

0 

. 8005 

4.  796 

-3.  545x10 

2 

.7995 

4.  789 

-3.  518x10 

4 

. 7990 

4.  786 

-3.  504x10 

6 

.7985 

4.  784 

-3.490x10 

8 

. 7980 

4.  780 

-3.476x10 

10 

. 7970 

4.  774 

-3.  449x10 

0.  1 

0 

.9128 

2.  742 

-9.  443x10 

2 

.9108 

2.  727 

-9.  252x10 

4 

. 9088 

2.  713 

-9.  066x10 

6 

. 9068 

2.  698 

-8.  885x10 

8 

.9048 

2.  684 

-8.  710X10 

10 

. 9028 

2.  671 

-8.  539X10 

0.  5 

0 

.9595 

1.  924 

-8.  335x10 

2 

.9550 

1.  880 

-7.  803x10 

4 

. 9505 

1. 840 

-7.  330x10 

6 

. 9460 

1.  803 

-6.  907x10 

8 

. 9410 

1.767 

-6.  484x10 

10 

. 9353 

1. 729 

-6.  049x10 

1.0 

0 

. 9723 

1.665 

-2.  064x10 

2 

. 9653 

I.  587 

- 1. 825x10 

4 

.9588 

1.  528 

-1.649x10 

6 

. 9510 

1.471 

-1.477x10 

8 

. 9410 

1. 409 

- 1.  297x10 

10 

. 9283 

1. 346 

-1.1 17x10 

1 
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* 

of  £ , 

* 

»•* 

Table  5-9. 

The  values 

v and 

for  LiCl  at  a 

range  of  molalities 

and  various 

It. 

m 

h 

* 

£ 

❖ 

V 

❖ 

AF£ 

0. 0001 

0 

. 5203 

17.  95 

-3. 840x1 0^ 

2 

. 5203 

17.  95 

-3. 840x10* 

4 

. 5203 

17.  95 

-3. 840x10* 

6 

. 5203 

17.  95 

-3. 840x10 

8 

. 5200 

17.  95 

-3. 831x10* 

10 

. 5200 

17.95 

-3. 831x10* 

0.  001 

0 

.6580 

9.  017 

-1. 185x10^ 

2 

. 6580 

9.  017 

-1. 185x10? 

4 

. 6575 

9.  014 

-1. 180x10? 

6 

.6575 

9.  014 

-1. 180x10° 

8 

.6573 

9.  013 

-1.  178x10? 

10 

.6570 

9.011 

- 1. 176x10° 

0.  01 

0 

.7975 

4.  777 

-3.463x10^ 

2 

.7970 

4.774 

-3.449X10' 

4 

. 7965 

4.  771 

-3. 435x1 0^ 

6 

. 7958 

4.  767 

-3. 415x10^ 

8 

. 7950 

4.  762 

-3. 395x10^ 

10 

.7945 

4.759 

-3. 382x10 

0.  1 

0 

.9048 

2.683 

-8.  709xl0*j 
-8. 538x10° 
-8. 372x10° 

2 

. 9028 

2.670 

4 

. 9008 

2.656 

6 

. 8988 

2.643 

-8. 211x10° 
-8. 016x10° 
-7. 846x10° 

8 

. 8963 

2.6  27 

10 

. 8940 

2.614 

0.  5 

0 

. 9420 

1. 771 

-6.561x10^ 

2 

.9365 

1.734 

-6.  135x10^ 

4 

. 9303 

1.696 

-5. 704x 10g 

6 

. 9228 

1.655 

-5. 248x 1 0^ 

8 

.9148 

1.616 

-4. 820xl0g 

10 

.9043 

1. 571 

-4. 333x10 

1.0 

0 

. 9455 

1.427 

, 10 
-1. 363x10 

-1. 184x10 

-1. 003x10 

2 

. 9335 

1. 366 

4 

.9185 

1. 301 

6 

. 8988 

1. 234 

-8.  211xl0g 

8 

. 8710 

1.  162 

-6.  3 5 5 x 1 0? 

10 

. 8215 

1. 070 

-4. 191x10 
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6.  Discussion 

It  is  interesting  to  observe  that  f , is  larger  than  f , and 

M corr)  * ±(  ex  p) 

increases  with  h.  This  is  to  be  expected  when  the  hydration  correction  is 

properly  applied  because  the  effective  concentration  is  larger  for  f , . than 

r 7 ' 1 ±(cxp) 

^°r  ^±fcorr)  leading  to  greater  interactions  which  depress  f , The 

±(  ex  p ) 

corresponding  free  energies  decrease  in  absolute  value  with  increasing  h 
reflecting  the  greater  stability  f salting-in)  for  lower  values  of  h.  Note  that 
we  are  referring  to  the  configurational  free  energy  alone  and  not  to  the  total 
one,  that  is,  to  that  part  of  the  free  energy  that  corresponds  to  activity 
coefficients  rather  than  to  activities. 

The  lattice  theory  yields  free  energies  that  undergo  inflexions  when  the 
mean  activity  coefficients  reach  their  minimum  values  and  reflects,  therefore, 
the  salting-out  process.  The  salting-out  terms  arise  in  the  derivative  of  the 
configurational  free  energy  in  relation  to  m from  the  derivatives  of  v and  of 
the  configurational  entropy. 

v,  in  addition  to  the  effect  of  r , reflects  hard-core  effects  ( Pitzer , 1973), 
ionic  field-dielectric  gradient  interactions  (Bahe,  1972)  and  ion-cavity  inter- 
actions (Desnoyers  and  Conway,  1964).  The  relative  contributions  of  these 
terms  cannot  be  assessed  at  present  because  the  above  authors,  after  demon- 
strating the  theoretical  validity  of  the  effects,  treated  them  individually  and 
obtained  the  relevant  coefficients  by  curve-fitting.  Furthermore,  Pitzer  used 
a Debyc-Huckel  type  term  for  non-specific  interactions  and  Bahe  as  well  as 


Desnoyers  and  Conway  introduced  a fixed  partial  degree  of  order  rather  than  a 
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variable  one  as  in  our  case. 

* 

The  degree  of  order  increases  and  v decreases  with  m as  is  expected 
from  the  theory  but  we  have  no  fir st- principle  explanation  for  the  increase  in 
order  in  going  from  LiCl  to  KC1  at  a given  m.  It  corresponds,  however,  to  the 
salting-out  trend  observed,  for  example,  by  Bahe  (1972). 
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Appendix  1 


Alternative  derivations  of  the  equation  for  f 


a)  Equilibrium  model 

Let  us  represent  by  C-a  a cation  on  the  anion  sub-lattice  (the  reciprocal 
lattice)  and  so  on.  At  equilibrium 


C-c  + A-a  £ C-a  + A-c 


(Al-1) 


fw  (C-a)(  A-c) 
fr2  (C-c)(A-I) 


(A  1 - 2) 


where  f and  f are  the  activity  coefficients  of  ions  in  the  reciprocal  and  in  the 
w r 

right  lattice  and  K is  the  equilibrium  constant. 

When  an  ion  is  added  to  the  right  lattice  (C  in  c or  A in  a)  the  energy 
increment  is 


AE  = 

r 


2 „*  * 

(ze)  A Ml  v 

D r 
e o 


(Al-3) 


where  asterisks  are  used  because  the  system  is  assumed  at  equilibrium.  For 


a wrong  ion  (C  in  a or  A in  ~) 


q 


A]  -2 


(I€)2  A t‘  v' 

4E  , 

r Dr 

e o 


For  reaction  (A  1-1)  then 


AE  = 2(AE  - AE  ) 
r-»w  w r 


4(ze  ) A £ v 

Au 

D r 
e o 


where  AE  is  the  replacement  energy. 

r-*w 

Setting  AE  = AF  as  AS  = 0 at  equilibrium, 


, 2 *2 
1 w 

AE°  = -kT  In  K = -k  T In  

^ f/o-wV 


At  the  standard  state  there  is  no  interaction  energy  as  the  standard  state  has 

the  properties  of  an  infinitely  dilute  solution  and  also  f = f =1.  Thus 

w r 


(1  - w ) 
o 


and  w =0.5,  the  definition  of  the  random  state,  w represents  the  fraction 
o o 


of  ions  in  their  reciprocal  lattices. 


i 


A I.  3 


As  K = 1 at  any  concentration,  from  Eqn.  (Al- 


2) 


In  addition. 


f * 

r w 


w 1 - w 


AE  = RT  In  f and  AE  = RT  In  f 
r r w 


w 


and,  therefore, 


f 


2 * '• 

r 2(ze)  A Z v 

RT  In  = AE  - AE  M 

V r w D r 


= RT  In 


w 


e o 


] - w 


s*c  y- 

As  £ = ( 1 - 2w  ),  w"  = 0-5(1  - l 


and 


v * 

. w I - / 

In  = In  - c 


1 - w 1 Jr  l 


* 


Thus 


. I,/  Z"AWZ<Mt’’ 


* aft 

v 


i + -e 


D r RT 
e o 


(A  1 -8) 


(AI-9) 


(A]  - 10) 


( A 1 - J I > 


D r k T 
e o B 


m 


( A ] - ] 2) 


A1  -4 


which  yields  Eqn.  (4-16)  when  r is  replaced  by  its  value  shown  in  Eqn.  (2-8). 

o 


b)  Roltzman  law  model 


Here  we  use  the  13oltzman  law  in  the  form 


probability  = constant  x e 


-E/kT 


(A  1 - 1 3 ) 


w , the  fraction  of  wrong  ions  at  equilibrium  is  also  the  probability  of  their 
occurrence.  Thus, 


(-AE  + AE  )/kT 
w r 


1 - w 


(A  1 - 14) 


AE  - AE 
w r 


1 - w 


2 (ze)2  Am  t v 
°e  ro  kBT 


(A  1 - 1 5) 


which  is  identical  to  Eqn.  (Al-12). 

Thus,  Eqn.  (3-16)  is  a necessary  condition  for  equilibrium  as£  and  v 

must  yield  minimum  values  of  A Fff,  that  is,  must  be  £ and  v . It  is  not 
7 c min  mm 

sufficient  because  only  when  LF„,  . , = AF;  , equilibrium  is  achieved.  Eqn. 

f(min)  (exp) 

(Al-12),  on  the  other  hand,  is  necessary  as  well  as  sufficient  for  equilibrium. 
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